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SATELLITE  MOTION  IN  AN  AX1 -SYMMETRIC  GRAVITATIONAL  YIELD 
PART  2:  PERTURBATIONS  DUX  TO  AN  ARBITRARY  Jj 

by 

R.  H.  Gooding 


SUMMARY 

This  Report  continues  th<  presentation  of  the  untruncated  orbital  theory 
begun  in  Technical  Report  8806'J.  The  effects  of  the  general  zonal  harmonic, 

,  are  now  covered,  the  main  results  being  a  trio  of  formulae  for  perturbations 
in  the  spherical-polar  coordinates  introduced  in  the  previous  paper  The  formu¬ 
lae  are  only  first-order  in  -,  but,  in  conjunction  with  the  second-order 
results  for  J2  published  in  Part  l,  the  complete  *et  formulae  may  be  regar¬ 
ded  as  constituting  a  second-order  theory,,  the  Earth's  J2  being  much  larger 
than  for  l  >  2  . 

The  mean  elements  of  the  theory  are  defined  in  such  a  way  that,  for  each 
,  the  coordinate-perturbation  formulae  have  their  simplest  possible  form,  with 
no  occurrence  of  zero  denominators.  The  general  formulae  are  used  in  a  rederiva¬ 
tion  of  the  results  for  J3  ,  given  in  Part  1,  and  in  a  derivation  of  results  for 
J«  • 

Numerical  comparisons  with  reference  orbits  are  held  over  to  a  later  report 
(Part  3) . 
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1  INTRODUCTION 

This  Report  Is  the  second  of  an  Intended  trilogy  devoted  to  satellite 
motion  about  an  axi-symmetric  primary,;  i.e.  about  a  gravitating  solid  of 
revolution.  Thus  it  continues  the  exposition  of  Ref  1 ,  which  will  henceforth  be 
referred  to  as  'Part  1',  Part  1  brought  together  the  principles  of  an  approach 
to  orbit  modelling  in  which  lengthy  expressions  for  short-period  perturbations  in 
the  usual  osculating  elements  are  compressed  into  concise  expressions  for  pert¬ 
urbations  in  a  particular  set  of  spherical-polar  coordinates;  it  then  proceeded 
into  the  presentation  of  a  complete  second-order  theory  for  perturbations  due  to 
the  zonal  harmonic  Jg  »,  and  a  complete  first-order  theory  for  J3  .  When  the 
primary  body  is  the  Earth,  J3  ’(and  every  subsequent  J*  )  is  of  order  J22  P/  so 
Part  1  may  be  regarded  as  describing  (for  J2  an<1  J3  only?  a  complete  second- 
order  theory  for  Earth  satellites,  where  'first  order'  refers  to  effects  of 
relative  magnitude  10"  3  .  Though  Part  1  has  only  recently  been  published,  a 
resumS2  of  the  theory  had  been  given  much  earlier. 

Two  other  papers  are  relevant  to  the  maturation  of  the  trilogy:  a  reient 
one3  on  mean  elements  (as  used  in  Part  1)(/  with  particular  reference  to  the 
relation  between  mean  semi-major  axis  and  mean  mean  motion;  and  a  much  earlier 
(and  more  important)'  paper*4,  of  slailar  title  to  the  trilogy’s,  that  established 
formulae  for  secular  and  long-period  perturbations  dae  to  the  general  (so 
general,  in  fact,  that  i  could  be  negative,  the  formulae  then  being  applicable 
to  lunisolar  perturbations)'.  The  present  Report  effectively  combines  the  new 
approach  of  Part  1  with  the  general  principles  and  notation  of  Ref  4,  the  result 
being  a  complete  theory  for  the  zonal  harmonics;  secular  and  long-period 
perturbations  are  applied  to  mean  orbital  elements,  and  short-period 
perturbations  to  coordinates. 

Part  3  of  tne  trilogy  will  be  largely  devoted  to  the  way  in  which  the  mean 
elements  evolve  over  periods  of  time  longer  than  just  a  small  number  of  orbital 
revolutions.  This  topic,  which  was  given  limited  attention  in  Part  1,(  is 
entirel}  neglected  in  Part  2.  It  is  intended  that  Part  3  will  also  give  details 
of  the  Fortran  program(s)  written  to  evaluate  the  accuracy  of  the  overall 
approach,  using  harmonics  up  to  Ji,  .  (The  variations  in  the  mean  elements  are 
computed  by  a  technique  that  involves  a  numerical  component  of  an  otherwise 
analytical  model,  aspects  of  this  technique  were  described  in  the  papcr5  that 
originally  outlined  the  author's  philosophy  of  coupling  a  hybrid  computational 
procedure  *0  the  coordinate-perturbation  approach.) 
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Other  authors  have  published  first-order  formulae  for  satellite 
perturbatio.13  due  to  the  geopotential;  they  usually  address  the  subject  more 
generally  than  here,  by  covering  the  tesseral  harmonics  as  well  as  the  zonal 
harmonics.  The  first  entirely  general  results  were  derived  by  Groves6, ;  in  an 
analysis  of  formidable  complexity,  whilst  the  classic  eference  is  the  text-book 
of  Kaula?.  The  very  generality  of  the  formulae  a  Refs  6  and  7  makes  it 
difficult  to  write  down  expressions  for  individual  effects,  however,,  and  it  is 
not  even  easy  to  show  that  the  two  se-.s  of  formulae  are  formally  equivalent  (the 
full  first-order  expression  for  the  perturbation  in  mean  anomaly  is  omitted  in 
Ref  6,  and  the  supplementary  terms  are  only  added  as  an  afterthought  in  Ref  7). 

Much  of  the  difficulty  in  the  general  analysis  arises  from  the  need,,  when 
covering  the  tesseral  harmonics  as  well  as  the  zonal  harmonica,  to  allow  for  the 
rotation  of  the  primary.  The  uniformity  of  this  rotation  with  time  makes  it 
natural  to  work  with  M  (mean  anomaly),  rather  than  v  '(true  anomaly),,  as 
integration  variable,,  but  this  inevitably  leads  to  infinite  summations.  When  the 
analysis  is  restricted  to  the  zonal  harmonics,  however,,  use  of  v  (rather  than  M  5 
leads  to  expressions  that  are  free  of  Infinite  summation,  and  Zafiropoulos6  has 
recently  published  untruncated  formulae  for  the  first-order  perturbations  in  the 
orbital  elements  due  to  the  general  Jg,  .  The  formulae  of  Ref  8  are  much  more 
explicit  than  those  in  Refs  6  and  7,,  but  this  is  unfortunately  at  the  expense  of 
some  very  long  expressions  -  it  takes  more  than  five  pages  to  express  the  basic 
formulae,  and  even  then  the  supplementary  terms  of  the  perturbation  in  M  are 
again  absent.  Now  it  will  emerge  from  the  present  Report  that  the  formulae  of 
Zafiropoulos  can  be  expressed  much  more  concisely.  The  real  breakthrough  comes,, 
however,;  when  the  short-period  perturbations  in  elements  are  replaced  by 
perturbations  in  coordinates.  If  it  were  not  for  the  rotation  of  the  primary, 
this  procedure  could  be  immediately  extended  to  the  tesseral  harmonics*;  for 
orbits  of  sufficiently  low  eccentricity  there  i3  no  difficulty,  and  very  simple 
general  formulae  were  given  in  Refs  5  and  9,  having  originally  been  derived 
during  o  stud/10  of  Navstar/GPS. 

As  w.  ;h  Part  1,,  a  List  of  Symbols  is  appended  to  the  Report;  it  is  almost 
entirely  consistent  with  the  List  of  Part  1,,  the  few  e:ceptions  being  noted.  The 
meaning  of  every  new  symbol  is  fully  specified  in  the  text,  but  only  minimal 


*  Appendix  A,  which  is  in  the  nature  of  a  postscript,  outlines  what  is  involved 
in  the  extension  for  a  non-rotating  primary,,  and  a  separate  paper  is  planned 
for  later  publication. 
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explanation  is  given  for  those  carried  over  from  Part  J.  This  is  true,,  in 
particular,  for  standard  symbolism:  thus  we  note,  straight  away,,  that  the 
orbital  elements  used  are  a,  e,  1,,  ft,  u  and  M  ,,  an  arbitrary  one  of  which  is 
denoted  (generieally),  by  c  ;  also  M  -  <j  ♦  /  ,'  where  /  is  s.iorthar.d  for 
|n  dt  ,  the  integral  being  taken  from  epoch  to  current  time.  We  continue  to  make 
use  of  the  quasi -elements  ip  ,,  p  and  L  ,,  reallj  only  defined  at  the 
differential  level;  thus,,  dip  -  dw  ♦  c  dft  (where  c  -  cos  1  )',,  dp  -  do  ♦  q  dip 
(where  q2  ■  1  -  )'  and  dL  -  dM  ♦  q  dip  . 

As  explained  in  Part  1,,  each  osculating  element,,  c  ,,  may  be  regarded  as 
the  sum  of  a  mean  element,  c  and  a  short-period  perturbation,  6c  >,  so  that 

C  -  ?  ♦  6c  •  (D 

A  ’semi-mean1  element,,  c  ,,  13  also  needed  (see  section  3*2  of  Part  1),,  but  in 
Part  2  we  will  usually  ignore  the  distinction  between  c  and  c  •  The  effect  of 
this  neglect  is  that  we  do  not  distinguish  between  the  quantities  denoted  by 
6c  t,  6pc  and  6c  In  Part  1,,  normally  using  6C  here  in  the  sense  of  6pc  of 
Part  1;  towards  the  end  of  the  Report  '(in  deriving  the  perturbations  due  to 
Ji;  ,,  in  section  8.5),  we  remind  the  reader  of  the  additional  terms  (split  beWeen 
C  and  6c  as  explained  in  Part  1)]  that  are  needed  to  express  '(f  lr3t-order), 
perturbations  in  full.  Not  even  the  distinction  between  osculating  elements  and 
mean  elements  is  of  significance  ir.  evaluating  the  right-hand  sides  of  equations 
in  general,  since  second-order  perturbations  are  not  taken  into  account  in 
Part  2,  but  the  following  important  distinction  (on  left-hand  sides)  between  c 
and  c  is  worth  noting:  Lagrange’s  planetary  equation  for  c  constitutes  the 
starting  point  of  analysis  for  the  element  c  »,  whereas  a  formula  for  c  is  part 
of  the  goal  of  that  analysis. 

The  analysis  is  greatly  facilitated  by  using,  instead  of  <*>  and  u 
(argument  of  latitude),,  the  modified  quantities  u’  and  u*  ,,  where 

u*  -  to  -  •J’1''  and  u'  -  u  -  in  .  (2) 

To  avoid  any  confusion,  it  is  remarked  that  the  use  of  the  accent  (prime  sign) 
here  has  a  connotation  entirely  different  frcci  that  which  distinguishes  a 
(osculating  semi-major  axis)  from  a'  ;  the  latter  quantity  is  an  absolute 
constant  of  the  motion  '(under  zonal  harmonics  only),  which  (as  shown  in  Part  1) 
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constitutes  the  best  choice  for  mean  serai-major  axis  (cf),  to  whatever  order  the 
analysis  is  conducted.  With  u'  now  introduced,  it  is  convenient  to  define  here 
the  much-used  quantities  Cj  and  Sj  ;  thus* 

Cj  -  cos  (jv  ♦  ku’)l  and  sj  •  sin  (jv  ♦  ku')'  •  (3) 

When  there  is  no  ambiguity  in  regard  to  k  ,,  the  superfix  (but  never  the  suffix)' 
will  often  be  omitted.  (Warning:  Cj  and  Sj  ,  as  used  in  Part  1 ,,  identify  with 
-C|„g  and  here.)' 

As  the  primary  is  assumed  axi-syrametric,  we  3tart  from  the  potential 

function  u/r  ♦  I  Up  »  where  the  individual  terms  of  the  disturbing  function  are 
t  * 

given  (in  the  usual  notation)  by 


ui  -  -  £  Jt  (R/r)1  Pt< sin  8)  .  U), 

The  value  of  i  in  the  summation  is  normally  taken  to  run  from  2  upwards,,  since 
the  cases  i,  -  0  and  l  -  1  are  essentially  trivial,,  but  the  general  formulae 
to  be  developed  cover  the  case  l  -  1  without  difficulty;  both  •tri-rial1  cases 
are  instructive  and  are  Interpreted  in  section  8,  following  Ref  4.  If  the 
concept  of  an  axl-symmetric  primary  is  generalized  to  allow  for  mass  outside  the 
orbital  region,  as  well  as  inside,,  then  (H)  can  be  extended  to  cover  negative 
i  ,  as  in  Ref  *4;  the  only  change  needed  in  the  expression  for  U^  is  that 
is  replaced  by  P-j.;  •  Our  overall  requirement  is  to  integrate  the  planetary 
equations  for  the  general  U^  ,,  thereby  obtaining  the  first-order  contributions 
to  each  f  and  6c  »,  and  then  to  combine  the  6c  (for  the  six  elements)  into 
6r,  6b  and  6w  ,,  the  corresoonding  perturbations  in  the  spherical-polar 
coordinate  system  attached  to  the  mean  orbital  plane;  the  latter  is  specified  by 
I  and  ft  K  and  the  transformation  from  the  |(r,b,w)-coordinate  system  to  the 
usual  rectangular  equatorial  system  is  described  in  detail  in  Part  1. 

In  section  2  we  decompose  as 

ut  -  I  (5) 


*  Thi3  notation  leads  to  more  concise  expressions  than  if  the  trigonometric 
argument  was  Jv  ♦  ku’  ,  as  was  originally  planned;  the  disadvantage  is  that 
the  Kepler-constant  quantities  are  now  and  S_|<  ,  rather  than  Cq  and  Sq 
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where  0  £  k  £  8.  and  u(£  is  only  non-zero  when  k  has  the  sane  parity  as  l 
(or  as  -  1 ,,  in  the  extension  to  l<  0)’  .  The  decomposition  arises  as  8  ,, 

In  ,,  1+  effectively  replaced  by  i  ,,  and  this  involves  the  Introduction  of 
certain  'lilies  of  inclination  functions.  The  functions  a£(I)'  were  originally 
introduced  in  Ref  4  and  are  used  again,,  but  quantities  A^  ,,  proportional  to  the 
Aj£(ii  values,  are  actually  more  convenient.  Related  functions,  and  asscciated 
quantities,  are  also  introduced,  and  recurrence  relations  are  given.  These 
relations  (and  corresponding  relations  for  the  eccentricity  functions,  referred 
to  in  the  next  paragraph)  are  required  here  in  the  development  of  the  theory,  but 
chey  are  also  important  as  computational  aids  in  the  implement at io*  f  the 
theory.  The  uj£  can  be  treated  separately  in  all  the  analysis  up  to  the 
derivation  of  the  6r  and  the  $w  ,  but  there  i3  a  complication  in  the 
derivation  of  ib  ;  this  will  be  handled  by  the  introduction  of  another  index, 

«  (/  which  is  always  of  the  opposite  parity  to  l  (and  hence  k)‘. 

Following  the  elimination  of  8  ,  we  must  also  eliminate  r  ,  using  the 
basic  formula  of  the  ellipse 


£ 

r 


1  ♦  o  cos  v  , 


(6) 


before  the  planetary  equations  can  be  integrated.  This  involves  families  of 
eccentricity  functions,,  which  are  introduced  in  section  3.  The  functions  B|(e) 
were  originally  introduced  in  Ref  4,  but  the  quantities  By  ,,  proportional  to 
them,  are  actually  more  useful.  Recurrence  relations  are  given,,  and  these  are 
even  more  important  than  the  relations  for  the  inclination  functions.  r.t  is 
implicit  in  the  use  of  the  By  that  every  uj[  could  be  further  decomposed, 
into  ■  utJ  say,  but  we  prefer  not  to  do  this,  postponing  the  introduction  of 
the  By  until  u£  ,,  in  each  planetary  equation,  ha3  been  eliminated  in  favour 
of  an  explicit  expression,  thus  the  notation  is  not  needed. 


The  development  of  each  planetary  equation,  via  first  the  A^  and  then 
the  By  ,  is  the  topic  of  section  4.  As  already  remarked,,  we  avoid  infinite 
summations  by  retaining  v  as  an  argument  of  the  equation  (as  opposed  to 
eliminating  it  in  favour  uf  M  ) ,  and  indeed  we  make  it  the  integration  variable 
by  applying  the  relation 


•  n  q*3  (p/r)2 


(7) 
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Each  equation  now  expresses  d  /dv  ,  rather  than  ?  ,  as  a  (finite)'  sun  of  terms 

k  k 

in  v  .  Each  such  term  is  just  a  multiple  of  Cj  or  Sj  ,  as  it  turns  out,  so 

the  integration  of  the  equation  is  immediate.  Terms  with  k  ♦  J  f  0  lead  to  the 

6c  by  definition.  Terms  with  k  ♦  j  -  0  ,  on  the  other  hand,  are  effectively 

constant  over  the  short  term,  and  contribute  directly  to  c  ;  when  k  -  0  (and 
so  also  J  -  0  ),,  the  Integrated  contribution  is  a  secular  perturbation,  whilst 
the  terms  with  k  j*  0  contribute  to  the  long-period  perturbation.  The 
distinction  is  Important  for  Earth  satellites,  because  the  secular  variation  of 
w  (due  to  J2)  oust  be  allowed  for  in  integrating  the  perturbation,  but  the 
subject  was  dealt  with  in  Part  1  and  will  be  picked  up  again  in  Part  3;  there 
will  be  no  further  reference  In  Part  2  to  this  coupling  between  Jg  and  the 
other  . 

Formulae  for  the  (£  due  to  Uj  )  are  collected  in  section  5,  whilst 

the  reduction  of  the  appropriate  dc  to  formulae  for  dr  ,,  db  and  dw  is  the 
subject  of  the  next  two  (and  much  longer \  sections.  As  described  in  previous 
papers  there  is  an  Important  distinction  between  the  integrations  required  for 
the  types  of  term:  for  the  ,,  the  process  leads  to  definite  integrals 
(see  Parts  1  and  3).  necessarily  zero  if  taken  over  zero  time  from  epoch;  in 
fck- components  of  the  d?  ,  on  the  other  hand,  the  process  leads  to  epoch- 
independent  indefinite  integrals  that  (apart  from  the  complication  of  semi-mean 
elements)  satisfy  (1).  But  indefinite  Integrals  contain  arbitrary  constants,, 
where  a  'constant'  in  the  present  context  is  any  quantity  that  is  Independent  of 
the  fast-varying  v  ,  i.e.  would  be  a  true  constant  for  motion  in  a  fixed 
ellipse.  It  is  only  when  these  constants  are  all  assigned  that  (at  the  first- 
order  level)  the  mean  elements,,  c  ,,  are  fully  defined. 

It  has  been  noted  that  enormous  advantage  accrues  from  taking  3  to  be  the 
exact  quantity  a'  (defined  by  the  energy  integral,,  as  explained  in  Part  1),,  but 
there  are  no  immediately  compelling  reasons  for  associating  particular  constants 
with  any  of  the  other  five  elements.  We  therefore  base  our  choice  on  the 
philosophy  of  making  the  expressions  for  dr  ,,  6b  and  6w  as  simple  a3 
possible.  These  expressions,  which  constitute  the  most  important  results  of  the 
Report,,  are  presented  in  their  general  form  in  section  6;  each  of  the  three 
expressions  involves  a  summation  over'  the  index  J,  with  the  integration  constants 
for  tne  elements  (other  than  a  )  not  yet  taken  into  account. 
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In  the  general  formulae  Just  referred  to,,  certain  values  of  j  In  the 
summations  would  Involve  terms  of  zero  denominator,  and  it  is  by  the  elimination 
of  all  these  terms  that  the  integration  constants  (other  than  for  a  )  are 
chosen.  This  is  the  subject  matter  for  section  7,,  which  completes  the  entire 
analysis.  An  outline  of  the  material  in  this  section  is  as  follows.  First,, 
the  formula  for  the  mandatory  constant  In  5a  (for  each  u£  )  is  recorded, 
essentially  as  a  matter  of  completeness.  Second,  the  constants  are  derived  for 
5e  and  5M  that  validate  the  omission  of  the  terms  with  particular  J  that 
would  otherwise  arise  in  5r  .  Third,  the  constants  are  derived  for  6i  and  58 
that  do  the  same  thing  for  5b  .  Fourth,  special  terms  (with  particular  J  ) 
in  5w  ,,  that  could  not  be  included  in  section  6  because  t».'*y  are  Induced  by  the 
constants  In  5e  and  SM  ,  are  obtained.  Finally,  the  constant  in  5u  (for 
earh  Uj  )'  is  derived. 

Examples  of  the  general  formulae  of  section  6,  together  with  the  special 
terms  in  <w  derived  in  section  7,  are  given  in  section  8:  first,,  for  the 
trivial  c  >sea  l  •  0  and  l  •  1  ,,  the  interest  in  which  '.as  been  remarked;  then 
for  2.-2  and  l  •  3  >,  leading  ,(as  a  useful  overall  checx)  to  results  already 
known  from  Part  1j  finally,,  fo*  l  -  4  ,  leading  to  formulae  not  hitherto 
published. 

2  FUNCTIONS  OF  INCLINATION  REQUIRED  IN  EXPANDING  THE  POTENTIAL 

Following  Ref  we  expand  P^(sin  8)  ,  required  in  (U),,  via  the  addition 
theorem  for  zonal  harmonics  (or  Legendre  polynomials);  thus 


Pt(sin  6)  -  l  uk  ff-;  kjj  Pj[(0)  P*(c)  cos  ku*  . 


:<8> 


Here  u0  -  1  ,,  uk  •  2  if  k  >  0  ,  and  the  Legendre  function  p£  is  defined  by 


k  dkP,(c) 

p*(c)  .  3<_i_ 


O) 
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The  second  factor  (the  k'th  derivative)'  in  (9)’  is  a  polynomial  in  c  , 
which  (with  k  S  A)‘  does  not  vanish  when  c  -  1  ,,  its  value  then  being 
(Jt  +  k)l/{2k  k!  (ft  -  k) ! }  .  Hence  this  factor  may  be  normalized*',  in  a  certain 
useful  sense,  and  we  write 


dck  2k  kl  (ft  -  k)l  * 


(10) 


where  a{( i )'  is  a  pure  polynomial  in  s  (-  sin  1 )'  if  k  has  the  same  parity  as 
l  ,  but  has  an  additional  factor  c  if  k  and  l  are  of  opposite  parity;  in 
each  case  the  constant  term  in  the  polynomial  is  unity,,  by  the  normalization. 
Explicit  expressions  for  the  a£(1)  are  given  in  Table  1,  for  values  of  A  and 
k  up  to  6. 

We  can  now  rewrite  (8)  as 

l 

Pt(sin  8)  -  I  <*tk  sk  Aj(l);  cos  ku’  ,  (11) 

k-0 


where  the  constant,-  aik  ,,  is  given  by 

ail<  -  uk  ?\W  /(2*  k!)'  .  '02) 

A  different  constant,,  Cj  ,  vas  used  In  Ref  1,,  Incorporating  a  factor  associated 
with  the  eccentricity  functions  of  section  3;  it  Is  given  by  -2  (  k  )  a^k  ,, 
where  (k)  is  the  usual  binomial  coefficient,,  m  here  (and  throughout  the  paper; 
being  used  to  denote  a  general  integer,  with  negative  values  allowed;  when 
l  <  0  ,  p£  in  '(12)  must  be  replaced  by  P.t-1  ,,  so  that  aik  -  ■  but  the 

relation  of  ajk  to  the  of  Ref  4  is  unchanged.  It  is  clear,  from  the  last 
paragraph,,  that  P^(0)  (or  P^.^O))  vanishes  when  k  and  l  (or  -A  -  1  )  are 
of  opposite  parity,,  and  it  may  be  shown  that  when  the  parity  is  the  same, 


*  This  ’normalization,  which  nas  nothing  to  do  with  th*  standard  normalization  of 
the  spherical  harmonics  and  their  J-coefficients,  leads  directly  to  Ak(i)  -  1  ,, 
one  of  the  pair  of  starting  values  for  the  recurrence  relation  (21).  For  some 
purposes  a  different  normalization  is  preferable  such  that  A£(i)  i3  defined 
for  all  i  2  0  ,  and  A&(i)  -  1  ;  the  family  of  normalized  functions  can  then 
be  extended  in  a  unified  manner  when  the  orbital  theory  is  to  cover  the 
tes3eral  harmonics  (Append!  A). 
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Pj(0)' 


i-U-k? 


(-1) 


2*  {+<«. 


a  *  k) i _ 

♦  k)}!  (iU  -  k))' 


03) 


In  substituting  (11)'  into  (4)  it  is  of  great  benefit  to  introduce  a  new*  quantity, 
AHk  defined  &y 


Aik  "  <Hk  sk  Ai(i)'  •-  !(14) 

where  p  (i  •  a(1  -  e^))  is  the  seml-latus  rectum  (or  parameter )  of  the  orbit; 
the  equation  applies  when  l  <  0  ,  so  long  as  the  suffix  of  A  is  replaced  by 
-t  -  1  .  The  use  of  A^  permits  us  to  write  the  general  term  of  (5)',  using 
also  the  notation  of  (3).  as 

-  -  £  (p/r)**;  Aik  eg  .  (15)' 

It  will  be  noted  that,  whereas  Aj(i)'  is  defined  and  useful  regardless 
of  parity,,  '(and  hence  )  is  only  non-zero  when  k  and  i  are  of  the 
same  parity.  The  zeroes  come  from  ait<  ,  for  which  the  non-zero  values,,  up  to 
k  ■  i  •  6  ,,  are  given  by  the  like-parity  entries  of  Table  2.  (The  Table  has  been 
extended  back  to  l  *  -7  ,  to  illustrate  the  identity  of  with  ^  when 

fc  <  0  .)]  However,,  a  use  will  be  found  for  quantities  that  behave  in  the  opposite 
way  from  and  Aik  ,,  anticipating  which  we  define  (with  bold  letters  to  make 

the  distinction) 


’  u<  (l  ’  K  *  Pjm(0)/(2<  <l  U  (16) 

and 

AIk  ’  Jt  (R/,>)<'  °U  s<  *l<‘>  •-  (,7> 

where  k  ha3  been  replaced  by  <  to  signify  tnat  we  now  have  quantities  that 
are  non-zero  only  when  <  and  l  are  of  opposite  parity.  Half  of  Table  2 
(for  all  %  )  is  devoted  to  oiK  ,  since  these  quantities  can  be  Included  with 
a^k  on  a  chequer-board  basis. 
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We  will  require  derivatives  of  the  inclination  functions.  It  Is  evident 
from  OO)  that 


U  -  k)(l  +  k  ♦  l)'  _  .k+1,.x 
- 2(K  >  1) -  3  A*  (1)  s 


from  this  and  (14)]  it  follows  that  the  (partial)'  derivative  of  with  respect 

to  i  Is  given  by 

*ik  ■  h  Hk  **''  {k<:  I  ir-l) f  *ru>}  ■  ('5>  i 

where  f  -  s2  .  The  quantity  In  (curly)  brackets  is  the  Dj£(i)]  of  Ref  4.  We 
will  also  require,  finally,  the  particular  combinations  of  A^k  and  A£k 
denoted  by  Ajk  and  Ajk  ,,  and  given  by 

4  •  ks'’  4  *•’  *ik  >  <2°> 

the  s"1  and  c"1  factors  do  not  imply  singularities,,  as  they  must  always 
cancel  via  k  A^k  ana  Ajk  respectively. 

The  Aj[(i)  and  (and  hence  the  A^k  )'  may  be  computed  with  the  aid  of 

recurrence  relations.  A  fixed  k  was  stipulated  in  Ref  4  for  the  formula 

(1  ♦  k)'  a£( i )'  -  (2i  -  1)'  C  aJ.jU)  -  U  -  k  -  1)'  A*J_2<i )  ,  (21) 


(21)'  is  even  valid  for  l  -  k  +■  1  ,,  if  an  arbitrary  (but  finite)  Ak_ ^  ( i )'  is 
assumed.  However,  it  is  usually  more  useful  to  stipulate  a  fixed  l  ;  the 
required  formula  was  giver  by  Merson11,  being 


A#(i)  -  c  Aa*^(i)  -  ^  7- k  ~  - 

c  Ai  Sl)  4(k  ♦  V) (k 


—  r  k\*2a) ,, 


valid  for  i  -  2  t  k  S  0  with  the  starters  A^( i >  -  1  and  a|~\(I)  -  c  ; 

(22)  is  also  valid  for  k  -  i  -  1  ,,  with  an  arbitrary  (finite)  A|+1(i)  .  Either 
of  the  two  preceding  ’pure’  three-term  recurrence  relations,,  (21)  or  (22;,  can  be 
used  with  just  one  ’mixed’  such  relation  to  generate  all  the  relations  connecting 
the  A*(i)  ;  perhaps  the  simplest  mixed  relation  (with  neither  l  nor  k  fixed) 
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U  ♦  k)  c  Aj(i>;  «  (1  -  k)  A^(i)'  ♦  2k  A^”1  (i)  .  (23) 

For  the  we  have  the  relation,  for  proceeding  along  a  'fixed  diagonal* 

of  Table  2  (with  l  >  0  and  a  constant  value  of  £  -  k  )', 


l  *  k  -  l 

alk  *  uk.,  k  “*-i,k-1 


(2k) 


whilst  to  proceed  to  a  lower  diagonal  we  have 


uk  (k  ♦  1)' 

“ik  "  rn< — a£-l  ,k+1  • 


(25) 


These  relations  suffice  to  generate  all  the  oAk  from  ao,o  ■  1  •  Similar 
relations  permit  the  generation  of  all  the  from  aito  -  -1  ,  they  can  be 

dispensed  with,  howev^'',  since  it  follows  from  (12)',  ,03)  and  06)'  that 

V  “  - Vl,<  “  '  *LT "  *  (26) 


Though  it  is  the  A^k  '(and  Atic  )  that  wt  actually  require  to  carry  through 
the  paper,  recurrence  relations  are  not  offered  for  these.  To  preserve  parity  if 
one  suffix  is  fixed,  it  would  be  necessary  to  use  alternate  values  of  the  other; 
there  seems  little  point  in  doing  this,  though  a  valid  relation  could  easily  be 
obtained,  for  example  by  applying  (21)'  three  times.  There  are  simple  relations 
between  the  Aik  and  A^  ,  however.  We  will  need  two  of  these  in  section  7.3, 
namely, 


and 


l 


(  At,k»1 

l  Vi 


A£,k-1 

Vi  . 


(27) 


l 


A£.k»l 

.  Vi 


A£,k-1 

Vi 


2A*k 
uk  ’ 


(28) 


Postscript,  in  regard  to  equations  (27) -(3D,  it  should  have  been  noted  that 
Aek/uk  and  Aik/Uk  actually  reduce  to  -fcc”1  Af  ,^/u,^  and  £c-1  A^  k+1'uk*1 
respectively,  results  that  are  implicit  in  the  a’nalysis  of  section  6.5,  the  uk 
factors  could  be  avoided  by  allowing  negative  k  and  *.  (see  the  footnote  of 
page  HQ), 
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these  being  true  for  1  s  k  *  l  (with  i  and  k  of  the  same  parity)'.  For 
k  -  0  we  only  have  one  relation,  given  by  direct  addition  of  (27)  and  (28)  '(and 
really  only  one,  involving  dirsct  subtraction,  when  k  «  Z  );  thus, 

*A*,1  "  ‘  *Ai,0  •  <29) 

We  could  use  (27)  and  (28)  to  obtain  expressions  for  Ajk  ,  defined  by 
(20),,  but  instead  derive  them  directly.  On  substituting  for  Aik  from  (14)'  and 
for  A^k  from  ,09)',,  then  in  forming  A j»k  we  find  that  the  term  in  a£(D 
cancels  out  and  we  get  ,(for  0  $  k  $  z  ,,  and  z  and  k  of  like  parity) 

Ajk  -  (B/p)»  atk  ■  »•’  ^  *r1(1>  •  (3°>' 

For  Ajk  ,  on  the  other  hand,  the  combination  of  aJ(1)  and  aJ+\(1)'  is  such 
that  {22),  with  k  replaced  by  k  -  1 ,:  is  immediately  applicable,  leading  to 
(but  only  for  1  $  k  ^  z  now)' 

Ajk  -  2k  Ja(R/p)*  aJlk  c-1  sX-1  aJ-1(1)  ;  {3D 

for  k  •  0  ,  this  would  give  a  false  value  of  zero,,  the  correct  result  being  the 
same  as  is  given  by  (30)',  since  a£  q  -  -A £  0  .  The  formulae  (30)’  and  (3D  are 
used  in  the  analysis  of  6b  in  section  6.2.  (See  also  the  footnote  to  page  15.)' 

3  FUNCTIONS  OF  ECCENTRICITY  USED  IN  THE  SUBSEQUENT  ANALYSIS 

The  term  U^  of  the  potential,  specified  by  (4),,  has  now  been  decomposed 
into  the  u£  defined  by  (15),,  the  latitude  (6)'  having  been  eliminated.  The 
longitude  was  absent  from  U^  from  the  beginning,  because  of  axial  symmetry,,  co 
it  remains  to  eliminate  the  radius  vector  (r),  .  This  can  be  done  by  appeal  to 
(6),  but  {as  noted  in  section  1)  we  will  in  practice  postpone  the  use  of  (6) 
until  the  setting  up  of  each  planetary  equation,  so  the  present  section  is 
preparatory  in  nature.  Further,  it  is  not  (p/r)*+1  ,  in  (15),,  that  must  be 
eliminated,  but  (p/r)*"1  ,  as  a  factor  (p/r)2  is  retained  to  effect  the  change 
of  integration  variaole  defined  by  (7D 

It  is  evident  from  (6)  that  an  expansion  of  the  form 
»- 1 

(p/r)*~l  -  l  uj  3m  cos  Jv  (32) 

J-0 
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is  possible,  for  12  1,  and  we  regard  B^j  as  defined  by  this  expansion; 
clearly, ,  B£j  is  a  polynomial  in  e  .  We  shall  find  it  useful ,,  and  entirely 
natural,  to  extend  the  definition  of  B^j  to  negative  j  ,  by  defining 
B*-J  "  Bjt,  | J  |  and  t0  take  Blj  *  0  when  |j|  *  i  .  On  this  basis,,  and  using  the 
notation  of  (3),  we  can  replace  ,(32)  by 

'(p/r)1"'  -  l  By  Cj  ,,  (33) 

where  the  summation  effectively  runs  from  J  -  -•  to  J  »  ,,  so  that  there  is 

no  need  for  explicit  summation  limits. 

To  make  use  of  some  results  from  Ref  4,  we  first  demonstrate  that  the  Bjj 
are  directly  related  to  the  Hansen  X  functions  of  classical  celestial 

mechanics^,  such  that 

By  -  q24*1  .  ,(3“1 

Hansen’s  functions  {of  eccentricity)]  are  defined  (uniquely)]  by  the  existence  of 
the  expansion,,  for  all  integral  l  and  J  ,,  regardless  of  sign,, 

(r/a)*  exp(\Jv),  -  l  exp(  unM)  ,,  (35) 

m  n 

where  •  -1  and  the  summation  runs  frc*n  -•  to  ♦“  .  Only  when  m  -  0  ,, 
which  is  the  case  with  which  we  are  concerned,  is  X*J  a  simple  (finitely 
expressible  in  elementary  functions)  function  of  e  (and  it  is  precisely  because 
of  this  that  we  change  the  variable  from  t  to  v  in  the  planetary  equations,, 
thus  avoiding  infinite  expansions  in  M  ) . 

To  demonstrate  (34),  we  first  replace  the  index  t  by  -(t  ♦  l;  in  (35), 
and  then  Integrate  over  a  revolution  of  M  .  We  get  (from  the  real  part  of  the 
result)' 


2* 

|  (p/r)l+1  cos  Jv  dM  -  2tj  q2(i+,)'  X5i_1’j  .  (36) 

0 
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But  from  (33)., 


2*  2* 

J  (p/r)i+1  cos  Jv  dM  -  |  l  Bjjj,  (p/r)2  cos  mv  cos  Jv  <JM  .  (37)’ 

o  o  a 

We  apply  (7)  to  change  the  integration  variable  to  v  on  the  right-hand  side  of 
(37);  only  If  m  •  ±j  do  we  retain  a  non-zero  term,  and  In  fact 

2it 

|  (p/r)*+1  cos  jv  dM  -  2ir  q3  .  (38) 

0 

Then  (3*1)  Is  Immediate  from  (36)  and  (38). 

Some  comments  related  to  the  notation  are  worth  makl  '.g  before  we  proceed 
further.  In  principle  we  are  reserving  the  suffix  k  for  the  A  functions  and 
J  for  the  B  functions,,  but  ii  section  5,  where  only  the  value  -k  arises  for 
J  ,,  we  will  naturally  encounter'  B4k  .  We  would  also  rather  naturally  change  the 
notation  from  J  to  k  In  (35)  if  we  w»re  following  the  traditional  path7  In 
which  the  Integration  variable  is  M  and  the  expansion  of  (15)'  is  by  (35)' 
directly.  After  replacement  of  l  by  ~U  ♦  1 }  ,,  the  Hansen  function  would  then 
appear  as  ,  which  is  nowadays  (following  Kaula7)  usually  expressed  (when 

1  5  0  )  as  Gj,pq(«)  ;  here  p  -  -fU  -  k)'  which  must  be  Integral  (assuming  l 
and  k  to  be  of  the  same  parity)',,  and  q  -  ra  -  k  .  Introducing  also  the 
notation  G^  ,  which  the  present  author?  has  recommended  as  preferable  to 
Gipq  •  we  raay  extend  (3*0  by  writing 

q-(2i-l )  By  .  x5t-c.j  .  -  oi,-j(e)  .  .(39) 

Gooding  and  King-Hele^  nave  recently  reported  o"'  the  G  functions  that  are 
relevant  to  resonant  satellite  orbits,  Ref  13  includes  the  listing  of  a  Fortran 
program  (by  Alfred  Odell)  that  computes  the  functions  for  arbitrary  values  of 
l  K  k  and  Kaula1  s  q  ,,  by  quadrature. 

We  can  now  use  the  identity  (3*0  to  tie  into  the  analysis  of  Ref  *4 .  Thus, 
we  may  express  the  e-polynomial  B^  K  when  2.  2  1  and  0  £  j  <  i  ,  in  terms  of 
a  normal i7ed  such  polynomial,,  the  connecting  relation  being 

BIJ  *  (*-  J  1 )  <e/2)J  B|(e)  ,  (HO) 
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B|(e)  in  x(#0)  is  a  polynomial  in  ,,  with  constant  term  unity  by  the 
normalization.  Explicit  expressions  for  the  s|(e)  are  given  in  Table  3,  for 
values  of  l  and  j  up  to  7  and  6  respectively.  There  is  an  evident 
resemblance  between  the  B-j[(e)  and  the  A^(i)  ,  a  significant  difference  being 
that  the  new  functions  run  from  t  -  1  and  not  l  -  0  ;  the  resemblance  is  not 
fortuitous,,  3ince  it  can  be  shown  that11 

Bi(8)  •  Ja(-  t  « li-itQ-1; .. 

from  which  it  follows  that 

b{(o)  -  ql-J"1  AjJ_,(taiT' le)  .  U2) 

In  contradistinction  with  the  A^,  ,,  however,  It  la  usually  much  better  to 
work  with  the  By  directly  (In  recurrence  relations,,  for  example),,  rather  than 
through  B{(e)  and  (40)'.  One  reason  for  this  la  that  only  alternate  values  of 
the  Alk  are  non-zero,,  whereas  (for  |j|  <  i  and  e  •  0)  all  the  By  are  non¬ 
zero.  Further,,  no  difficulty  arises  with  the  By  when  J  <  0  (as  already 
noted,  and  see  also  Appendix  B),  whereas  b{(b)  would  then  be  Infinite  .< If 
| j  |  <  j.) .  «e  can  even  allow  l  to  be  negative  (or  zero)  as  well  as  J  .  The 
validity  of  this  follows  from  the  universality  of  (34)  -  the  universality  Is 
brought  out  by  Table  4,  which  lists  By  for  t  running  from  -3  to  +4  and  J 
from  -1  to  *3« 

The  entries  In  Table  4  form  triangular  blocks  of  Tour  types.  First,  Tor 
l  >  0  and  |j|  <  l  ,,  we  have  the  quantities  that  can  be  expressed  by  (40)  when 
J  10.  Secondly,  for  l  >  0  and  |j  |  i  l  ,  we  have  (two  blocks  of)  zeroes. 
Thirdly,  for  ISO  and  |j|  S  -i  ,  we  have  quantities  that,,  when  J  1  0  ,  can  be 
expressed  by  a  formula  complementary  to  (40),,  viz 

BtJ  -  (*  J  1)  (e/2)J  q2*-l  B^,(e)  ;  (*3) 

a  formula  equivalent  to  this  was  given  in  Ref  4,  the  application  being  ,(as  noted 
in  section  1)  to  secular  and  long-period  perturbations  associated  with  exterior 
(rather  than  interior)  mass.  Finally,  for  l  £  0  and  |j|  >  i  »  we  have  .(two 
blocks  of)  quantities  that  are  most  conveniently  expressed  in  terms  of  8  (not 
now  denoting  latitude,  a3  previously)  and  q  ,  rather  than  e  and  q  ,,  where 
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*  -  TT7  ’  (M* 

no  attention  wis  paid  to  these  quantities  (or  their  equivalents),  in  Ref  4,  there 
being  no  application  for  them,  but  the  formula  for  Bqj  is  derived  here,  for 
completeness,  in  Appendix  B.  (Other  entries  in  this  last  pair  of  blocks  are  then 
derivable  frrm  recurrence  relations.)'  Before  leaving  Table  4,  we  note  that  the 
formula  for  the  X  or  G  function  corresponding  to  Bjy  is  immediate  from  the 
Table,  in  view  of  (34);  tnus  it  is  only  necessary  to  apply  the  factor  q1"2*  , 
which  introduces  a  negative  power  of  v,  when  there  is  not  one  already  present 
and  cancels  it  out  when  there  is! 

In  regard  to  derivatives  of  the  eccentricity  functions,  it  can  be  shown  (by 
working  from  (41)),  and  easily  verified  from  Table  3  that,,  for  1  S  j  <  i,  , 

•^■1b|(o)|  •  2Je-’  .  («5) 

The  universal  formula  for  the  derivative  of  is 

Bij  •  <»  -  '>  et-u-i  -  j  »■'  •  <«> 

For  1  i  J  <  l  ,  this  follows  from  (45);  for  general  entries  in  Table  3.  it  can 
be  verified  with  the  aid  of  q’  and  3*  ,  which  may  be  expressed  as  -e/q  and 
S/eq  respectively.  However,,  because  we  only  introduce  the  By  after  each 
planetary  equation  has  been  set  up,  we  effectively  only  use  (46)  in  expressing 
the  rates  of  change  of  tne  mean  elements.  Since  this  involves 

■k (q  Aix  Btk>'  ■  »■’  ik  'q2  bJk  ♦<«-'> *  Btk>  •  cn 

we  define 

Elk  "  S2  Bik  *  -  O  «  ;  (k8) 

then  (46),,  re-expressed  via  the  :  ecurrence  relation  (56),  leads  to 

E£k  -  e"1  (JLe2  -  k)  B4k  ♦  (t  -  k)  .  (49) 
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(By  symmetry,  there  is  a  parallel  expression  for  E^k  that  involves  Blk  and 
b£  k+1  •)  Table  5  gives  explicit  expressions  for  the  E^k  ,  with  t  running 
fromtr3  to  +4  as  in  Table  4;  only  the  entries  in  which  k  has  the  same  parity 
as  1  (or  -A'-  1  if  i.  <  0  )  are  useful  in  practice,  and  entries  for  k  >  t  (or 
-1+1  if  %  2  0  )  are  omitted  entirely  (for  k  i  l  >  0  they  would  all  be  zero). 
The  Ejik  are  related  to  the  Ej[{e)'  of  Ref  4  by 

Etk  "  6-1  <e/2)k  (*  jj  M  E^e>  *.  (50), 


when  i,  >  0  ;  for  1  2  0  ,,  the  extra  factor  q2*-4-1  is  required  (cf  (43)',  where 
the  additional  factor,  in  relation  to  (40),,  is  q2*-*1  )'. 

The  By  may  be  computed  from  recurrence  relations  for  the 
will  now  be  given,  but  in  developing  the  theory  it  is  more  useful 
relations  for  the  By  themselves,  so  these  will  also  be  given. 

(  i  0),  the  recurrence  formula  (from  Ref  4)  is 

U  ♦  J  -  1)'  Bj[(e),  -  (21  -  3)  8|.1  -  u  -  J  -  2)  q2  Bj.2(e)'  ,  (50 

valid  for  1  l  J  +  3  with  the  starters  Bj+1(e)’  and  both  unity.  For 

fixed  t  (  l  3).  on  the  other  hand,  the  formula  is11 


B{(e)  ,,  which 
to  have  such 
For  fixed  J 


Bi<e> 


It  •  .)  ■  2)(1  <  J  t  I) 
MJ  *  U(J  ♦  2) 


a2  B{‘2(e)'  , 


(52), 


valid  for  l  -  3  2  J  2  0  with  the  starters  B^_1(e)'  and  B^'2(e)  again  both 
unity.  The  resemblance  of  (51)  and  (52)  to  (21)'  and  (22)  respectively  follows 
from  the  remark  leading  up  to  (42). 

The  recurrence  relations  for  the  By  ,,  that  correspond  to  (51)'  and  (52)., 
respectively,  and  are  valid  for  all  l  and  j  ,  are  (when  symmetrically 
expressed) 


lit  -  !>'  q2  Bt.,  j  -  1(2 1  -  O'  *  U2  -  J2)  Blt1|J  -  0  (53) 

and 

<j  -  1)  e  Byj.,  ♦  2J  By  ♦  (J  ♦  1)  e  ByJ+1  -  0  ;  (54) 
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the  mandatory  symmetry  in  (5*0 »,  to  satisfy  the  unchanging  value  of  By  under 
the  operation  j  ♦  -j  ,  is  evident. 

As  remarked  for  the  inclination  functions,  either  of  the  ’pure’  relations,, 
(53)  and  (54),  can  be  used  with  a  'mixed’  relation*  to  generate  all  the 
recurrence  relations  connecting  the  By  .  Here  the  mixed  relations  are,  in 
particular,  those  that  connect  three  out  of  four  of  the  By  lying  'around  a 
square'  of  index  duplets;  if  the  square  consists  of  the  duplets  (1,  J),y 
(i  '•  1,  J),  (i,  J  +  1)  and  (t  ♦  1,  J  ♦  1),  then  the  four  mixed  relations 
connecting  them  (all  of  which  we  shall  require  in  the  sequel)  are 

t  By  -  U  -  j)  ByM  *  la  BtiJ+1  -  0  .  (55J 

tq2  By  -  (l  -  J),  Byyj  ♦  ( l  *  j  ♦  1  )'  B  Byyjy  -  0  ,  (56) 

*«  0lj  *  1  -  U  *  J  ‘  D  Byyjy  -  0  .(57) 

and 

U  -  J)  e  BtM>J  ♦  lq2  Bl>J„  -  it  *  J  *  1)  Bu1>Jy  -  0  .  (58) 

If  we  re-order  the  terms  in  the  last  two  relations  and  replace  j  by  J  -  1  ,,  we 
get  relations  which  are  symmetric  pairings  of  (55)  and  (56)',  viz 

l  By  -  (t  ♦  J)'  Byt#J  ♦  i  e  B^j.j  -  0  (59)' 

and 

lq2  By  -  U  ♦  J)'  By,tJ  ♦  ( t  -  J  ♦  1  ?  <!  Byyj.,  -  0  .  ( 60) 

Of  this  set  of  relations,  (55)  ar.d  (59)]  can  be  obtained  at  once  from  (54)  and  the 
relation  equivalent  to  one  given  (for  the  Hansen  functions)]  by  Zaf iropoulos®,,  viz 

te  -  h,y i?  -  Bfi,j  ; 

this  is  of  a  different  ’shape’  from  our  triangles-around-the-square  relations, 
but  is  perhaps  the  simplest  recurrence  relation  of  all. 


*  Note  added  in  proof:  Ref  19  indicates  that,  for  inclination  functions,  pure 
relations  are  computationally  preferable  to  mixed  relations  (see  also  Ref  20). 
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4  RATES  OF  CHANGE  OF  OSCULATING  ELEMENTS 

Ir.  tnis  section  we  use  Lagrange’s  planetary  equations  to  develop  the  rate 
of  change  of  each  of  the  orbital  elements  (a,  e,  i,,  0,  w  and  M)  due  to  ,  the 
term  of  the  disturbing  function  specified  by  (15).  Each  rate  of  change  is  to  be 
with  respect  to  v  ,  rather  than  t  ,,  expressed  as  a  finite  trigonometric  series 
(assuming  %  >  0  ,,  as  we  now  always  do,  except  in  section  8.1),,  with  v  as  the 
variable.  The  v-independent  terras  of  each  dc/dv  are  then  isolated;  they 
effectively  contribute  to  the  time  rate  of  change,  t  »,  of  the  mean  element,, 

X  expressions  for  the  c  being  held  over  to  section  5-  The  remaining  terms  of 
dt/dv  can  at  once  be  integrated  to  provide  contributions  to  the  short-period 
perturbation,  5r  .  The  result  of  the  integration  is,  in  fact,;  so  ’immediate' 
(apart  from  the  question  of  the  integration  'constants'  already  referred  to  in 
section  1)  that  we  will  not  bother  to  write  down  formal  expressions  for  the  five 
it  other  than  5a  this  Is  to  emphasize  the  fact  that  it  is  the  combinations 
of  the  5c  into  5r,  5b  and  5w  that  are  of  interest  (being  the  topic  of  section 
6),,  not  the  5C  themselves. 

The  perturbation  5a  is  a  special  case  because  it  can  be  obtained  without 
integration.  As  in  Part  1,  however,,  we  also  derive  da/dv  from  the  appropriate 
planetary  equation,  as  a  prototype  for  the  derivation  of  the  other  dc/dv  .  By 
bringing  in  the  quasi-elements,,  $  and  p  ,  it  is  possible  to  develop  each 
equation  in  terms  of  the  partial  derivative  of  u£  with  respect  to  a  single 
quantity. 

H.1  Seoi-raajor  axis 

As  in  Part  1,,  there  is  an  absolute  constant  of  the  motion,  which  we  denote 
by  a'  ,  such  that 


a  -  a’  (1  ♦  2aU/p)  ;  (62) 

this  is  an  exact  relationship  for  any  time- Independent  disturbing  function,  U  , 
and  in  particular  for  the  axi-symmetric  .  It  follows  that  there  is  no  long¬ 
term  variation  in  a  ,,  to  whatever  order  of  magnitude  the  perturbation  analysis 
is  conducted.  Further,  the  short-period  perturbation,  5a  ,,  is  given  exactly,,  on 
substituting  for  ujf  from  (15),  thus 

Sa  •  -  2a' (l'2  AJk  (p/r)1*’  .  (63) 
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This  does  not  mean  that  an  exact  perturbation  can  be  written  down  for  semi-major 
axis,  however,  as  the  right-hand  side  of  (63)  is  expressed  in  terms  of  osculating 
element;",  as  soon  as  mean  elements  are  introduced,  the  result  is  no  more  than  a 
first-order  perturbation  expression,  as  with  any  other  5  . 

To  present  6a  in  the  form  appropriate  for  use  in  section  6.1 K  we  combine 
Cq  with  one  of  the  factors  p/r  .  Thus, 

6a  -  -  aq"2  Alk  (p/r 5*  («c51  ♦  2Cq  ♦  eC*)  .  (64) 

We  retain  another  p/r  factor  explicitly,  and  expand  the  remaining  (p/r)*"1  by 
(33).  By  this  means  the  term  2Cq  .  for  example,  in  (64)  is  effectively 
transformed,  for  each  J  K  into  Cj  ♦  C_j  .  But  each  pair  of  terms  (such  as 
this)'  for  positive  J  ,  in  the  infinite  summation  of  (33).  is  matched  by  the  same 
pair  (in  reverse  order)  for  negative  J  ,  so  we  can  express  the  result  of  the 
expansion  as 

6a  -  -  aq"2  A^k  (p/r)]  l  Bjj  (eCj.,  ♦  2Cj  ♦  eCj+1)  .  (65)' 

We  now  develop  an  expression  for  da/dv  ab  initio,  using  the  general 
procedure  that  involves  the  planetary  equation  for  a  .  This  equation  is 


da  _  2  31) 

dt  *  na  3M 


and  on  substituting  for 


ujj  we  get 


da 

dt 


-  2naq"2  Aik  Cq}  ■ 


(66) 


(67) 


The  M-dlfferentiatlon  is  immediate,  since  p/r  is  given  by  (6)'  and  3v/3M  is 
q"3  (p/r)2  (cf  (7)).  We  transform  from  dt  to  dv  (again  using  (7)),  and  all 
this  leads  to 

—  -  aq-2  *4k  (p/r)1  |(k  -  l  -  1)  e  S_,  ♦  2k  S0  ♦  <k  ♦  i  *  I )  e  S, 1  .  (68) 
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Finally,  we  split  (p/r)5-  into  (p/r)5-"1  and  p/r  ;  then,  applying  (33)  and 
taking  advantage  of  the  matching  of  terms  for  positive  and  negative  j  ,,  we  get 

-  iaq'2  I  By  l(k  -  1  -  1)  e2  Sj_2  ♦  2(2k  -  l  -  1),  e  Sj_, 

♦  2k(2  ♦  e2)  Sj  ♦  2(2k  *  I  *  H  i  *  (k  *  l  *  1)'  e2  Sj,2)  .  (69)' 

This  may  be  regarded  as  a  prototype  for  all  the  dc/dv  ;  in  addition, 

(69)  is  used  in  the  derivation  of  de/dv  in  section  4.2.  The  equivalence  of 
this  result  with  the  v-derivative  of  6a  obtained  by  the  special  procedure  may 
be  verified,  most  easily  by  the  v-differentiation  of  ,(64). 

In  dealing  with  the  subsequent  c  we  will  be  isolating  the  component  of 
dc/dv  that  leads  to  secular  and  long-period  perturbations.  We  know  that  for 
da/dv  this  component  must  be  2ero,  both  from  the  special  procedure  and  from  the 
form  of  (66)  (since  a  term  of  U  that  is  free  of  short-period  variation  must 
tautologously  have  zero  M-derlvatlve)',,  but  it  is  instructive  (as  part  of  the 
prototype  for  the  other  c  $  to  obtain  this  result  from  (69).  The  terras 
independent  of  v  in  the  overall  J-sura  are  the  terms  in  S-k  '(  »  S^)  .  Since 

*  8i,J,  the  combination  of  all  such  terms  involves  the  factor 

(k  -  l  -  1)  e2  Btik.2  «  2(2k  -  t  -  1)  e  *  2k(2  ♦  e2)  Bu 

♦  2<2k  ♦  l  *  1)  e  *  (k  *  l  *  1)  e2  Bt  K.2> 

and  it  follows  from  three  applications  of  (54),  that  this  is  zero. 

4.2  Eccentricity 

We  develope  the  perturbation  in  e  by  first  obtaining  the  perturbation  in 
p  ly  since 


p  -  q^  a  -  2ae  e 

and  the  planetary  equation  for  p  is  Just 

dp  _  2q  3U 
dt  *  na  3u 


(70) 


(7D 
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On  substituting  for  u£  we  get 

P  -  -  2naq-1  Alk  (p/r)4*'  (C0)'  .  (72) 

Transformation  of  the  integration  variable  to  v  ,,  by  (7).,  yields 

-  2kp  Alk  (p/r)4"1  S0  ;  :<73> 

hence  we  get,,  from  (33)'  and  the  usual  argument  concerning  positive  and  negative 

J 


$  -  21<!>  *ik  l  SJ  •  <W 

We  get  the  long-term  variation  by  setting  j  -  -k  ;  thus 

P  -  2knp  Aik  Blk  S.k  .  (75) 

The  expression  for  S  is  now  immediate  from  $70)',  since  $  -  0  ,  but  it  is  not 
given  here,  as  the  complete  list  of  the  ;  is  given  in  section  5. 

For  the  v-derivative  of  the  short-period  perturbation,  $e  ,,  we  have  all 
the  terms  with  J  ♦  k  *  0  in  the  expression  given  by  the  combination  of  (6y)  and 
(73).,  according  to  (70)'.  This  combination  leads  to 

37  -  »  »tk  l  BtJ  (<k  -  t  -  1)e  Sj_2  *  2<2k  -  l  -  1)Sj_, 

*  6ks  Sj  *  2(2k  *  2  »  OSj.,  *  (k  *  l  ♦  l)e  Sj»2)  .  (76) 

As  already  indicated,  we  will  not  write  down  the  expression  for  6e  ,, 
involving  Cy2  Qtc».  8iven  immediately  by  integration  of  '(76) ;  immediate,,  that 
is,,  apart  from  the  ’constant  term’,;  in  C.k  f  tnat  we  are  not  yet  in  a  position 
to  assign.  This  term  effectively  replaces  the  infinite  term  (in  C.k  )'  that 
would  arise  if  we  had  not  removed  the  S.k  term  from  (76)  in  advance. 
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4.3  Inclination 

We  can  develop  the  perturbation  in  i  from  the  perturbations  in  p  and 
pc2,,  since 

d(pc2)/dt  -  c2p  -  2pcs  i  (77) 

and  the  planetary  equation  for  pc2  is  Just 

«<p£>.  .  iSS-i"  (78, 

dt  na  3P  * 

But  uj  is  independent  of  longitude,,  and  hence  of  Q  ,,  so  pc2  is  an  invariant. 
Using  (77) ♦,  therefore,  we  have  t  at  once  from  ,(75),;  whilst  41  will  be  based 
on  the  expression  for  dl/dv  derived  from  (74)',,  viz 


f  -  kcs-T  A*  I  By  .  (79) 

(Since  contains  sk  as  a  factor,,  there  will  never  be  a  non-zero  multiple 

of  an  uncancelled  s’1  .)' 

4.4  Right  ascension  of  the  node 

The  perturbation  in  (5  ,oraes  from  the  planetary  equation 


do  ^  1  3U 

dt  na2qs  3i 


On  substituting  for  u£  we  get 


fl  *  -  nq*3  s”’  (p/r)**1  Ajk  Cq 


(80)' 


'(80; 


We  low  apply  ,(?)  and  (33)  as  usual,,  getting  an  isolated  contribution  to 
ft  ,  together  with  the  expression  for  dQ/dv  ,,  viz 


d a 
dv 


3"1  Ajk  l  By  Cj  . 


(82) 
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4.5  Argument  of  perigee 


Introduction  of  gives  us  a  one-term  planetary  equation,  since 


d£  .  Q  5U 
dt  na2e  3e 


:C83), 


The  e-derlvatlve  is  much  the  most  complicated  of  the  partial  derivatives  of 
u£  ,y  since  e  is  an  argument  of  each  of  the  four  factors  on  the  right-hand  side 
of  (15).  Thus  we  get 


i>  -  -  ne"1  Q'2  (p/r),1*1  Coj, .  '( 84) 

But 

.  ai.iu-h,,  (ss) 

and,,  using  the  expressions  for  3r/3e  and  3v/3e  (equations  (41)]  and  (42)  of 
Part  1),, 

-^■|(p/r),**1  Coj  -  q'2  (p/r)**1  (U  ♦  1)(cos  v  -  a  -  e  sin2v)]  Co 

-  k  sin  v  <2  *  e  cos  v)  Sq}  (86) 

so  that  (84)'  reduces  to 

i  -  ne^q'3  Aik  (p/r)**'  {k  sin  v  (2  ♦  e  cos  v)  So  - 

U  ♦  1)]  cos  v  (1  ♦  e  cos  v)’  Co}  .  (87) 

We  make  the  standard  expansions  of  the  trigonometric  products  In  (87),,  and 
then  apply  (7)  and  (33)  as  usual.  This  leads  to 

-  -it'1  *tk  l  B(j  its.  »  1  -  k)e  Cj.2  ♦  2(1  ♦  1  -  2k)Cj_, 

*  2(1  ♦  1)e  Cj  *  2(t  *  I  ♦  2k)Cj,,  *  (S  *  1  »  k)e  Cj.2)  .  (88) 
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To  get  ft  ,,  we  pick  out  the  coefficient  of  C_k  .  Thus 


i  •  -  i  ne-'  Alk  (U  *  1  -  k)e  8lik_2  ♦  2(t  *  1  - 

♦  2 (t  *  Ort  Blk  *  2U  *  1  ♦  2k)Blik„  ♦  U  ♦  1  ♦  k)e  Bt>k+21  C.k  .  (89) 

But  this  can  be  simplified  by  three  applications  of  (5^)'»,  which  lead  to 


$  -  -i.  "2  Alk  (<te2  -  k)8lk  ♦  <t  -  k)e  Blik.,l  C.k  .  (90) 

We  can  now  introduce  the  quantity  E«k  ,  to  get  a  concise  expression,,  since  by 
(1*9)’  we  have 

ft  -  “  ne-1  Alk  Eik  C_k  .  (91)' 


To  get  &  and  the  appropriate  terms  of  dw/dv  ,  we  combine  (91)'  and  (the 
residual  terms  of)'  (88)'  with  ft  and  (82), ,  respectively,  using 

w  -  4'  “  c&  .  (92) 

H.6  Mean  anomaly 

We  start  by  studying  p  ,,  since  our  final  cne-term  planetary  equation  is 


d£  _2__3U 
dt  ”  na  9a 


(93) 


Remembering  that  A^k  ,  in  ,,  is  itself  a  function  of  3emi-maJor  axis,  we 
obtain 

i>  -  -  2U  *  t)nq'2  (p/r)1*'  Alk  C0 

and  hence 

It  ■  - 2(1  *  ,)!i  stk  l  •  :<95> 

In  particular, 

ft  -  -  2U  ♦  1  )nq  Aak  0tk  C,k  ,  (96) 


T«  89022 


30 


and  from  (91)'  we  therefore  also  have 


i  •  -  nq  Atk  U(t  ♦  1)Btk  -  e"1  Elk|  C.k  ,, 


since 


o  -  p  -  . 


From  (88)  K  similarly,;  the  v-derivative  of  5o  Is  given  by  the  v-dependent 
of 


$f  -  *  e"'<I  Atk  I  By  (U  ♦  1  -  k)e  Cj.2  *  2U  ♦  1  -  2k)Cj_, 

-  6(t  ♦  1 )e  Cj  ♦  2(1  ♦  1  *  akJCj*,  *  ,U  *  t  *  k)e  Cj,2}  . 

But 

M  •  c  ♦  /  , 

where  (with  t  standing  for  time) 

t 

/  -  |  n  dT 
0 

and  (assuming  only  Uik  to  be  operating)’ 

n  -  a'  -  3nq’2  Atk  Cp/r )*■*'  C0 
by  (63)  and  Kepler's  third  law. 

From  $101 )  and  (102)  It  follows  that 

—  .  y  .  3,  Atk  l  By  Cj  ,, 
by  the  u3ual  procedure.  We  nay  then  write 

t  -  n’  *  3nq  A,k  Blk  C.k  . 


(97) 


(98) 

terms 


(99) 


(100) 


(101) 


(102) 


(103) 


(104) 
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from  which  ft  is  available  on  combining  with  (97).  Finally ,,  on  combining  the 
residual  terms  (those  with  k  +  J  *  0  $  of  (103)’  with  (99),  we  find  that  the 
v-derivatlve  of  6M  is  given  by  the  v-dependent  terms  of 

—  •  *  e'’q  Atk  l  BtJ  |(t  +  1  -  k)e  Cy?  *  2(1  ♦  1  -  SkJCj., 

-  6(t  -  1  )e  Cj  ♦  2tt  ♦  1  ♦  2k)Cj,t  *  ll  *  1  ♦  k)e  CJ4,21 

To  conclude,,  we  note  that  a  very  much  simpler  result  than  (105)  is 
available  for  the  non-singular  5L  .  Thus  from  (88)  and  (105)  we  get 

w  •  - (2t  -  ’>«  s  ■ 

5  SECULAR  AND  LONG-PERIOD  ELEMENT  RATES 

In  this  section  we  collect  the  expressions  for  the  rates  of  change  of  the 
mean  elements,,  i.e.  the  ;  associated  with  .  As  we  have  seen  in  section  4, 
this  simply  amounts  to  listing  the  components  of  the  dc/dv  that  are  multiples 
of  either  S*k  or  .  When  k  »  0  ,  the  rate  of  change  is  secular;  for 
k  >  0  ,  it  is  long-period.  We  shall  not  be  concerned  with  the  build-up  of  actuaj. 
perturbations  from  the  ;  since  this  is  fully  dealt  with  in  Parts  1  and  3; 
suffice  it  to  say  that  there  is  no  difficulty  in  the  secular  perturbations,  but 
that  (even  in  a  first-order  analysis)  difficulties  arise  with  the  long-period 
perturbations,  in  particular  due  to  the  singularities  associated  with  zero  e  and 
zero  s  . 


(105) 


(106) 


Another  point  must  be  mentioned  before  we  list  the  c  •  As  the  expressions 
arise  from  terms  in  dc/dv  ,  but  were  treated  (in  section  4)  as  if  from  terms  in 
d£/dM  ,,  each  $  produces  a  hort-period  component  of  the  perturbation  in  c  ; 
i.e.  a  contribution  to  5;  is  induced.  These  contributions  may  be  amalgamated 
into  components  of  6r  ,  4b  and  $w  ,,  as  done  for  J3  in  Part  1  (section  7)  ■ 

The  issue  relates  to  the  definition  of  semi-mean  elements  (section  3.2  ibid)’,, 
which  is  outside  the  scope  of  Part  2;  it  should  be  clear,  frem  equations  (120)  - 
s(122)  m  section  6,  however,  that  no  difficulty  arises  in  the  amalgamating. 
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In  the  list  of  the  c  that  follows,,  we  note  that  the  maximum  value  of  k 
is  l  -  2  ,  since  8^  -  -  0  .  We  attach  an  explicit  subscript  (Ik)  to  each 
t  ;  then  our  first  result  is 


*4k  -  0  .  (107) 

For  e  ,,  it  follows  from  (70)  and  (75)!  that 

4  -  -  *"«'V  *tk  B»k  4  •  <10S> 

There  will  always  be  a  positive  power  of  e  to  cancel  the  factor  e*l  ,  it  will 
be  noted,  coming  from  k  Bik  . 

For  i  K  similarly,,  it  follows  from  (75)  and  (77)  that 

tlk  ■  kncs-1  Aik  Bik  s\  .  (109) 

Here  there  will  always  be  a  positive  power  of  s  ,,  coming  from  kAik  ,  to  cancel 
the  factor  s”1  . 

For  n  ,,  our  analysis  of  (8D  gives 

5  fyk  ■  "  n  Aik  Bik  c-k  ’  <no> 

The  formula  is  expressed  in  this  way,,  with  a  factor  s  on  the  left-hand  side,  to 
avoid  the  possibility  of  an  uncancellable  s_1  on  the  right-hand  side.  For 
long-period  r  vturbations,  there  is  a  singularity  difficulty  here,,  which  can  be 
dealt  with  33  indicated  in  Part  1  (section  3.5).  For  secular  perturbations  (and 
here  is  our  first  r on-zero  Z  when  k  -  0  )  there  is  no  problem,  since  in  the 
expression  for  Ajk  ,  given  by  (19),,  A*(i)  appears  with  the  multiplying  factor 
k  ,,  and  Aj[*\(i)  with  the  factor  f  . 

For  w  ,  we  use  the  final  result,  ( 91 ) ,,  for  $  .  Then  from  (92)  and 
(110)  it  follows  that 

«>Elk  -  R(«c  AJk  Blk  -  ,  Alk  Elk>  .  (ini 
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Again  the  fornula  is  expressed  like  this  to  make  the  right-hand  side  non¬ 
singular;  and  again  (because  Eak  contains  a  factor  e  when  k  -  0  K  as  seen 
from  Table  5)  there  is  no  difficulty  with  secular  perturbations. 

For  M  ,  we  combine  the  results  for  ?  and  7  -  n*  ,,  given  by  (97),  and 
(104)'  respectively;  thus 

e  -  nq  Aj,k  (Eik  -  (21  -  l)e  Blk}  clfk  .  (112) 

From  (48),  this  may  also  be  written  as 

e  -  nq3  Aik  Bjk  C*k  .  (113) 

From  the  definition  of  L  ,,  we  may  also  combine  (112)  with  (91);  this  gives  the 
non-singular  result 


-  -  (2*  -  1)'  nq  Alk  8lk  C*k  .  (114) 

As  usual,  a  factor  e  can  be  cancelled  froa  both  sides  of  (113)  when 
k  -  0  .  However,  there  is  a  simpler  way  of  dealing  with  secular  perturbations  in 
M  ,,  as  indicated  in  Part  1;  Ref  3  was  largely  devoted  to  this  topic,  and  the 
rest  of  this  section  conforms  with  the  account  therein. 

The  basic  idea  is  that  we  represent  the  secular  perturbations  in  mean 
anomaly  by  modifying  the  value  of  the  mean  mean  motion.  In  view  of  (113),  in 
fact ,,  we  write 


TV  -  n' ( 1  ♦  e'V  l  ,  BJ#0)  .  (115) 

where  the  summation  is  now  on  i.  ,  and  we  have  set  Oq  to  unity.  This  agrees 
with  equation  (11)'  of  Ref  3,  since  A^  Q  here  may  be  identified  with 
-JjCj(R/p)*Aj(i)  from  that  paper. 

The  logic  for  using  a’  a?  mean  semi-major  axis  (S)  is,  as  we  have  seen, 
compelling,  so  if  h  •  n*  ,  we  do  not  retain  Kepler's  third  law  in  its  simplest 
form.  This  is  of  no  consequence,  however,,  and  we  simply  write,  from  (115), 

n2s3  •  u(l  *  2e"'q3  l  A10  B'l  0)  .  (116) 
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Ref  3  gives,  as  equation  (15),.  an  explicit  version  of  '( 1 1 6 )  with  (even) 
values  of  l  up  to  8;  the  two  versions  can  easily  be  verified  as  equivalent  if 
we  note  that  {in  Ref  3)  is  Just  a  normalized  form  of  B[  Q  such  that 

Bi,0  *  ‘  1)()l  *  2>efy  •  (117), 

The  following  recurrence  relation  was  given  in  Ref  3: 

<1  -  DQj  -  (21  -  -  a  -  4>;  q2  Qt.2  j  i(  118); 

this  is  valid  for  £  2  4  K  with  Qg  ■  0  and  Q3  -  1  .  The  relation  may  be 
obtained  from  (117),  and  (46)',;  together  with 

(l  -  1)(i  -  3)  -  (1  -  2)(2t  -  5)  Bi-2,l 

-  u  -  2)(t  -  3<  <t2  Bt.3>1  ;<119)' 

which  derives  from  (53)  on  replacing  (£,  J)  by  (t  -  2,  1). 

It  is  very  convenient  that  Qg  *  0  .  It  means  that  for  first-order 
analysis  associated  with  Jg  (the  dominant  harmonic  of  the  geopotential)'  TT  is 
the  same  as  n'  . 

6  PERTURBATIONS  (SHORT-PERIOD)  IN  COORDINATES  -  GENERAL  CASE 

In  this  section  we  develop  general  expressions  for  the  5r  ,  5b  and  6w 
that  can  be  derived  from  the  first-order  5;  via  the  formulae  (taken  from 
section  3-3  of  Part  1) 

5r  •  (r/a)  5a  -  (a  cos  v)'  5e  ♦  (aeq'l  sin  v)  5M  ,,  (120) 

5b  •  (cos  u*)  51  ♦  (s  sin  u')  58  (121) 

and 

5v  -  6*  ♦  (q~2  sin  v  (1  ♦  p/r)}  5e  ♦  q~3  (p/r)2  5H  .  022) 

Special  cases  (derived  from  the  choice  of  integration  constants  in  the  5C  )]  are 

reserved  to  section  7,  but  in  counting  the  number  of  terms  associated  with  the 
general  we  have  regard  to  the  basis  on  which  these  constants  are  chosen. 
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The  d;  are  available  at  once  as  the  v-integrals  of  the  expressions  for 
d£/dv  in  section  4.  Generation  of  the  expressions  for  dr  and  dw  is 
essentially  straightforward  in  that  the  analysis  starts  with  the  dc  due  to  uj[ 
(which  d?  we  can  denote  by  6;^  )'  and  finishes  with  dr^k  and  dw^k  .  with 
db  ,  however,  there  is  a  complication,  due  to  the  appearance  of  u’  in  (121)',  as 
opposed  to  v  in  (120)'  and  (122);  as  already  noted  in  section  1,,  we  deal  with 
the  difficulty  by  deriving  dbt<  ,  rather  than  dbj,k  ,,  where  <  has  values  of 
opposite  parity  to  those  of  k  . 

We  do  not  give  expressions  for  dr  f<  db  and  dw  K  but  (as  is  clear  from 
Part  1)'  these  are  immediately  available  from  tne  expressions  for  dr  ,  db  and 
dw  ,,  Just  by  replacing  Sj  and  Cj  by  (respectively)  (k  ♦  j)'  n  Cj  and 
-(k  ♦  J)  n  Sj  .  We  can  do  better  than  this  if  we  allow  for  the  (overall)  rate  of 
change  of  <3  ,’  replacing  (k  ♦  J)  n  by  (k  ♦  J)'  n  +  k  w  ,,  assuming  Cj  and  Sj 
still  to  be  shorthand  for  Cj  and  sj  . 

6.1  The  perturbation  dr 

We  have  to  apply  (120)'  with  da  ,  de  and  dM  given  by  (65)  and  (the 
integrals  of)'  (76)'  and  (105).  We  find  that  the  integrals  combine  in  a  very 
natural  way,,  as  a  result  of  which  we  can  write  (with  dr  short  for  drlk  )' 

5r  -  (r/aj  6a  -  t  a  »tk  !  By  {.  ♦  3  k  ’'jcj.,  II 


,  2k  -  l  -  1 

1  .  2k  ♦  l  ♦  1 

,  k  ♦  l  -  t  ,  k  ♦  t  *  1  L  \ 

■l  X  ♦  J  •  1 

l  k  ♦  J  ♦  1 

■J  Cj  *  e 

3  k  *  J  k  *  J  *  2jCJ*') 

The  simplest  way  to  Incorporate  (65)'  is  to  note  that  this  can  be  decomposed 

into 


S  «a  -  -  -t  a  [  Bjj  (9(5-*  *  3 


*2[2HfH--2H-j4r)cJ  •  :<124>  II 

By  this  trick,  we  can  combine  (123)  and  '( 124)  at  once,  to  get,,  say,, 

«r  -  -  ;  a  |  By  Rj  .  <125) 
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It  can  be  seen  that  (125)]  is  a  summation  in  which  Rj  ,,  as  given  by  ](126)\ 
has  three  components;  each  component  is  expressed  as  the  sum  of  two  multiples  of 
the  same  'C  quantity'.  Let  us  separate  the  first  multiple  from  the  second  (in 
each  component  of  Rj  )#  feeding  them  back  separately  into  the  summation  of 
(125),,  so  that  we  have  two  distinct  summations  that  we  can  denote  by  I.  and 
£+  .  Thus  £_  involves  l  Rj_  ,  where 


k  *  J  -  2 


♦*¥rr^f  ej  *  3‘ 


Now  we  have  seen  (in  section  3)  that  all  sums  over  can  be  regarded  as 

running  from  -•  to  .  It  follows  that  we  can  rearrange  the  three  sets  of 
terms  in  IB^j  Rj_  such  that  (with  J  now  used  in  a  different  way)' 

I  BtJ  Rj.  •  l  (k  ♦  )  -  U-’Kj  *  t)e  *  2( 2J  *  l  -  t)Btj 

*  3(J  -  t>«  Bt>j.,|Cj  .  (128) 

We  now  invoke  the  recurrence  relation  (51*)  to  eliminate  •  so  that 

the  quantity  in  curly  brackets  in  (128)  becomes 

2<J  ♦  l  -  1)8tj  ♦  2(J  -  !)e  Bjj.,  ,, 

and  then  simplify  further,  using  (58)  with  both  l  and  J  reduced  by  1 ,  to 
reduce  this  to  2(1  -  1)q2  .  (We  get  the  same  result  oy  using  (5*0  to 

eliminate  first,,  and  then  simplifying  further  via  (56).)  Thus 

I  BtJ  Rj.  -  2(i  -  1)q2  l  (k  *  J  -  I)'1  Cj  .  '(129) 
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Similarly, 

I  Bjj  Rj+  -  -  2(i  -  t)q2  l  (k  *  i  *  I)'1  Bt.,  j  Cj  .  0  30) 

The  final  result  we  require  now  follows  from  (125),,  (129)  and  (130). 

Because  of  its  importance,  we  write  Cj  in  full.  Thus 

5r*k  -  -  U  -  1$  p  %  j  (j<  ;  3 -  7  3  .  ^  B^., tJ  ct.3  (ku'  *  Jv)  .  031)' 

Equation  (131)]  provides  a  general  formula  for  6r  due  to  U*  ,,  valid  for 

til.  (This  restriction  on  1  has  been  operative  from  the  beginning  of 
section  4.)  In  view  of  the  fact  that  k  only  takes  non-negative  values  of  the 

same  parity  as  1  K  it  should  be  noted  that  J  takes  all  values,,  but  with 

Bt-l,j  only  non“2ero  |j|  i  1  -  2  .  (This  applies  if  12  2;  but  the  case 

1  -  1  is  trivial  because  there  is  an  overall  factor  1  -  1  . )] 

If  J  -  -  k  t  1  ,,  there  is  a  zero  denominator  in  (130'.,  and  terms  with 
these  values  of  J  must  be  excluded  from  the  formula;  they  are  associated  with 
the  terms  in  4;/dv  that  were  hived  off  In  the  generation  of  the  t  •  In 
section  7  we  shall  determine  constants  for  <5e^k  and  such  that  the  terms 

with  these  two  values  of  J  are  forced  to  zero.  It  will  be  noted  that  all  the 

cosine  terms  occurring  in  030',  for  a  given  and  all  possible  k  ,  are 

distinct,,  except  that  if  k  ■  0  (l  even)'  then  equal  and  opposite  values  of  J 
lead  to  identical  terms  in  cos  Jv  . 

We  use  the  remarks  in  the  last  paragraph  to  provide  a  pair  of  formulae  for 
,  the  total  number  of  terms  required  to  express  4r  for  a  given  value  of 
1  .  One  formula  applies  when  1  is  odd,,  the  other  when  l  is  even.  In  both 
cases  the  number  of  J  values  for  each  k  (regardless  of  the  excluded  values,  if 
any)  is  21  -  3  if  12  2. 

If  1  is  odd,  there  are  $(l  ♦  1)  possible  values  of  k  ,  so  a  priori  the 
value  of  Nir  is  ±(t  ♦  0(21  -  3)  if  12  3*  But  this  must  be  reduced  by 
the  number  of  excluded  values  of  the  duplet  (k,  J)  .  If  k  -  1  ,  J  cannot  be 
-k  i  1  ,  so  there  is  no  value  to  exclude.  If  k  «  1  -  2  J  car,  (a  priori)  be 

-k  ♦  1  and  this  value  must  be  excluded.  If  k  £  1  -  4  ,  it  will  always  be 
necessary  to  exclude  both  -k  *  1  and  -k  -  1  .  Thus  the  total  number  of 
exclusions  is  1  -  2  .  Subtracting  this  from  the  a-priori  value,  we  get 
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*tP  -  *2  -  j(3l  -  1)  .  (132)' 

Values  for  l  up  to  15  (including  N1  r  -  0  ,,  which  i3  the  correct  value,  even 
though  the  above  analysis  only  applies  for  l  >  3  )  are  given  in  Table  6. 

If  l  is  even,  there  are  j-l  ♦  1  possible  values  of  k  ,,  so  a  priori  the 
value  of  N^r  is  ♦  2)(2l  -  3).  The  exclusions  are  as  recorded  before, 

amounting  to  i  •  1  now  if  l  s  4  but  it  is  also  natural*  for  N£r  not  to 

count  the  'duplications'  that  arise  when  k  -  0  ;  there  are  %  -  3  of  these 

duplications  if  t  £  H  ,  viz  for  2  <  |j|  <  i  -  2  (we  cannot  'discount'  for 

| J j  -  1  ,  since  both  values  have  already  been  'excluded').  Thus  the  total  number 
of  exclusions  is  effectively  increased  to  21  -  H  ,,  and  this  is  the  right  number 
even  when  1-2  (not  covered  by  the  argument  that  applies  for  1  >  4  only). 
Subtracting  this  value  from  the  a-priorl  value,  we  get 

Ntr  -  1 2  -  |(3t  -  2)  .  :<133>, 

Table  6  gives  values  for  1  up  to  16.  It  is  remarkable  that,,  as  a  result  of  the 
discounting  of  the  duplications,  we  have  a  formula  that  is  so  close  to  what  the 
improper  use  of  (132),  would  give,  the  value  by  ,0  33)'  being  larger  by  Just  j-. 
Further,  if  we  did  not  discount,  the  formula  for  1  >  4  ,,  viz  -  £(1  ♦  4)  , 
would  give  1,  instead  of  the  correct  2,  when  1  -  2  . 

It  is  noted,  in  conclusion,  that,  due  to  the  multiplier  (r/a)  of  6a  in 
(120),  it  would  not  be  a  simple  matter  to  null  the  'constant'  terms  of  dr  with 
a  choice  of  1  *  a’  ,  but  that  in  any  case  we  would  prefer  not  to  make  such  a 
choice.  Further,  the  constants  in  6a  and  6r  are  not  the  same,  partly  due  to 
the  multiplier  of  6a  referred  to,  but  mainly  to  the  way  in  which  the  terms  in 
de  and  6M  combine.  For  even  l  ,,  we  will  have,  in  particular,  a  coefficient 
of  Cg  (-  1)  equal  to  'U  -  DpA^  q  .  (See  '(159)'  for  the  constant  in  6a.)' 

6.2  The  perturbation  db 

We  get  db  from  '(121),  where  dl  and  dQ  are  given  by  the 
Integrals  of  179)  and  (82).  This  is  on  the  assumption  that  6b  {•  db^)  is 
associated  with  ujf  ,,  following  the  decomposition  of  by  (5).  We  shall 
shortly  find,  however,  that  Jt  is  much  more  convenient  to  decompose  the  total 


*  In  software  in  particular,  we  would  rather  double  a  computed  quantity  than  have 
to  compute  it  again;  In  the  general  analytical  formula,  (130,  however,,  there 
is  no  easy  way  to  indicate  a  special  situation  when  k  -  0  and  j  *  0  . 
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6b  (associated  with  )  as  £  50^  ,  where  the  summation  is  for  values  of  < 
that  are  of  opposite  parity  to  t  and  we  no  longer  associate  the  individual  6b 
(-  )  with  specific  components  of  . 

Xn  relation  to  ,  we  get 

-  - 1  h3  { Fry  cj  *  Fry  %  sj  3l"  “’}  ■  '(1W  I 

The  trigonometrical  products  can  be  replaced  by  sums,,  in  the  usual  way,  and  we 
can  then  invoke  the  notation  of  (20)  to  write 

Jbtk  -  -  ic  l  By  (k  ♦  JJ-'UJk  cj*1  ♦  AJk  cf'>  .  (135) 

This  expression  may  be  contrasted  with  (125)  and  (126),  for  6r  .  In  view  of  the 
difference  in  superfix,  as  well  as  suffix  (which  alone  varied  in  the  terms  of 
Rj  )’,,  in  the  two  C  terms  of  ,(135),,  we  would  now  like  to  combine  a  pair  of  terms 
with  different  k  indices,,  before  the  summation  over  the  J  index  operates.  We 
note  that  Ajk  and  ,  though  under  the  summation  sign  in  (135)',  are  actually 
independent  of  J  . 

With  the  philosophy  Just  referred  to,  we  make  the  new  decomposition 

«bt  ■  l  5btK  •:  S' 36) 

where  each  ib  '(•  $b^k  )  la  of  the  form 

Jb  -  I  Tj  BtJ  Cj  (137) 

and  we  require  an  expression  for  Tj  (or  Tilcj  to  display  all  the  index 
parameters).  We  note  first  that  since  (for  non-trivlal  results)  k  runs  from  0 
or  1  to  i,  (taking  alternate  values),,  it  follows  that,,  in  principle,,  <  runs 
from  -1  or  0  to  i  ♦  1  (again  alternate  values,  but  of  opposite  parity  to  k): 
for  the  minimum  value  of  <  ,,  only  the  term  in  Ajk  ,  in  (135),;  contributes  to 
Tj  ,  whilst  for  the  maximum  value  of  <  ,,  only  the  term  in  Aj^  contributes;  for 
intermediate  values  (if  any),,  both  terms  contribute.  But  we  can  straight  away 
dismiss  the  ’maximum  value’  (k  -  t  ♦  l),  because  A^  is  just  a  multiple  of 
s*  ;  from  this  it  follows  that  A£t  ,  defined  by  (20),,  is  zero.  (Also  -  0 
anyway!)  We  shall  find  that  we  do  not  require  the  'minimum  value’  (<  -  -1) 
either . 
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To  evaluate  Tj  In  general,,  we  use  (30)  and  '(31)’  for  A£k  and  Ajk  ,, 
respectively.  It  is  fortunate  that  we  require  the  first  with  k  -  <  -  1  and  the 
second  with  k  -  <  ♦  1  ,  since  this  means  that  we  pick  up  the  same  inclination 
function,,  Aj(i)  ,,  for  both;  moreover,,  it  is  aesthetically  satisfying  to  have  a 
direct  application  for  the  inclination  functions  for  which  subscript  and 
superscript  are  of  opposite  parity,  as  opposed  to  merely  an  application  in  the 
propagation  of  like-parity  functions*.  We  have 


-  h 


{ 


K  *  1 

K  ♦  J  ♦  1 


afc,«+1 


a  -  <  *  ou  ♦  <) 
n<{<  ♦  j-  i) 


s*  Aj(i)  .  036) 


The  quantity  in  curly  brackets  in  (138)  is  j,  pure  constant,,  in  which  the 
a.k  are  given  by  '(12) j  thus  the  first  a  involves  P|+\(0)  and  the  second 
involves  Pj-1(0),  ,,  these  being  given  by  (13).  By  relating  these  to  Pj+1(0)  ,, 
we  may  express  the  aforesaid  quantity  (after  some  algebraic  reduction)  as 


U  -  K  ♦  D‘ 
2<+1  <! 


i  (0) 


V-1 

.<♦3*1 


•c  +  J  -  1 


l 


But  Pjtl(0)’  is  related  to  at<  by  ( 1 6) ,  and  thus  to  by  (17).  Hence 

(138)  gives 


T  *  /'  U<»1  u<-1  \ 

J  ■  2u*  \<  ♦  J  ♦  1  "  <  ♦  J  -  1  J 


(139) 


The  preceding  is  ’general'  in  that  it  applies  for  1  £  <  S  l  -  1  '(or,,  more 

precisely,-  with  2  as  lower  bound  when  l  is  odd);  further,,  if  <  i  2  we  can 
obviously  cancel  the  three  appearances  of  u  .  We  still  have  to  cover  the  cases 
<  -  0  (Z  odd)  and  <  ■  -1  (i,  even)’,,  in  which  (in  principle)'  only  the  first  term 

in  curly  brackets  is  to  be  taken.  To  count  only  non-zero  terms  when  Tj  is 
substituted  in  (137)',,  the  restriction  on  J  is  that  |jJ  £  t  -  1  (cf  an  upper 
bound  of  1-2  in  the  analysis  for  r  );  we  shall  be  excluding  the  values 
J  -  x  t  1  ,  of  course. 


*  This  strengthens  the  view  (noted  in  other  papers,  and  in  section  8  of  Ref  13 
in  particular)  that  V  (of  either  parity)  is  a  much  better  index  parameter  than 
Kaula’s  p  (where  2p  -  Z  -  k  ,  referred  to  in  section  3  here).  When  the 
analysis  includes  the  tes3eral  harmonics  (Ref  9,  and  see  Appendix  A  here),  k 
takes  negative  values  (with  |k|  $  l  )  as  well  as  positive,  but  the  factor  uk 
in  (12)  is  not  required. 
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When  k  -  0  ,,  we  require  just  2/(J  ♦  1),  from  the  curly  brackets  in  039); 
but  for  each  J  >  0  ,  half  this  quantity  may  be  combined  with  half  the  corres¬ 
ponding  quantity  for  J  <  0  ,  to  give  1/(J  ♦  1),  ♦  1/(-J  ♦  l)  ,,  which  can  be 
rewritten  as  1/(j  ♦  1)  -  1  /(J  -  1)’.  It  follows  that  (139),,  with  the  three 
occurrences  of  u  deleted,  again  gives  the  right  results  (counted  separately  for 
j  >  0  and  J  <  0).  The  modified  formula  may  be  seen  to  apply,,  finally,,  for 
j  -  0  . 


When  <  «  -1  ,,  the  position  is  more  complicated.  First,  the  expression  by 
(139)  is  not  even  legitimate  now,,  since  k^K  is  not  defined  for  <  <  0  ;  the 
illegitimacy  arose  in  the  substitution  for  Aj^  since  (3D  does  not  apply  when 
k  -  0  .  But  (20)  indicates  that  a£  q  -  -AJ  0  ,  and  this  suggests  that  we  can 
relate  the  required  term,  involving  with  <  -  -1  ,,  to  the  te-m  in  ukC-1 

when  <  ■  1  .  Since  Cj*  -  CD  ,,  the  relating  will  involve  the  transposition  of 
positive  and  negative  values  of  J  ,,  and  this  is  also  necessary  to  identify 
u<4,i/(k*J*1)'  for  «  •  -1  with  uK_j/(<*J~1)'  for  x  •  1  .  In  short,,  we  can 
deal  with  <  -  -1  Just  by  doubling  the  second  term  in  curly  brackets  in  '(139) 
that  is  associated  with  <  •  1  .  This  means  that,,  yet  again*,  we  get  the  "ight 
result  from  '(139)  if  we  cancel  the  three  appearances  of  u  . 

We  can  now  write  down  the  final  result  we  require,;  on  substituting  (139)' 
into  (137)  and  expressing  Cj  in  full.  Thus 

5btK  -  -  t  *u  I -fTryrTTTTT-p-ry  8tj  oos  Uu’  *  Jv>'  •  !,"o) 

As  already  Indicated,,  this  formula  is  unlike  '(13D’,;  the  corresponding  one 
for  fir  ,  in  that  it  cannot  be  taken  in  isolation  as  relating  to  a  sub-component 
of  .  It  is  like  d3D]  in  one  respect,  however,,  in  that  terms  of  fib^  with 

J  ■  -k  t  1  are  excluded.  In  section  7  we  shall  determine  constants  for  fii^ 
and  fifl^  (k  ,,  not  <  ,  now  being  the  appropriate  symbol)'  such  that  these  terms 
are  forced  to  zero. 


*  The  universality  of  this  procedure  (cancelling  the  u  )  stems  from  the 
original  introduction  of  uK  into  the  definition  of  olK  .  if  we  dispensed 
with  this  factor,  but  used  positive  values  of  <  as  well  as  negative  ones 
(see  also  Appendix  A),  then  we  would  find  nothing  special  about  the  /alues 
tl  and  0  in  the  first  place. 
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We  proceed  to  obtain  a  pair  of  formulae  for  N£b  the  total  number  of 
terms  (without  duplication  of  Cj  )]  required  to  express  6b  for  a  given  value  of 
l  Whether  l  is  odd  or  even,  the  number  of  J  values  for  each  <  (not  dis¬ 
counting  excluded  values),  is  21  -  1  ,,  correct  for  all  £  U  1)  this  time. 

If  l  is  even  (which  we  have  seen  to  be  the  simpler  case),,  there  are 
possible  values  of  <  ,  so  a  priori  the  value  of  N£b  is  ±1(21  -  O'  .  When 
k  ■  i  -  1  j  can  be  -k  +  1  but  not  -<  -  i  ,,  so  there  is  Just  one  value  to 
exclude.  For  all  other  <  ,  values  of  •<  ±  1  are  both  possible,  so  the  total 
number  of  exclusions  is  l  -  1  .  Subtracting  this  from  the  a-priori  value,,  we 
get 

N£b  -  l2  -  i(3l  -  2),  .  0«U 

Interestingly,,  this  Is  the  same  as  N£r  given  by  (133)-  Values  for  l  up  to  16 
are  given  in  Table  6. 

If  l  is  odd,  there  are  ±U  ♦  D]  possible  values  of  <  so  a  priori  the 
value  of  is  ±(t  ♦  0(21  -  O'.  There  is  again  a  single  exclusion  if 

<  -  i,  -  1  P/  and  two  otherwise,-  so  there  are  i  basic  exclusions  (assuming 
1  i  3  ) .  in  addition,,  however,;  there  are  l  -  2  duplications  when  <  -  0  ,,  and 
these  can  be  discounted  for  N£b  (though  not  for  3b  itself  -  see  also  the 
footnote  in  section  6.1)',  so  the  effective  number  of  exclusions  is  21  -  2  \ 
this  value  applies  even  when  i.  -  1  .  Subtracting  this  total  number  of 
exclusions  from  the  a-priori  values  we  get 

Nlb  "  *2  "  *(a  ' 

which  is  one  more  than  for  the  corresponding  N£r  .  Values  for  l  up  to  15  are 
given  in  Table  6. 

In  conclusion,  it  is  worth  remarking  that  if  the  planetary  equations  are 
used  in  Gauss's  form,,  as  opposed  to  Lagrange's,  {and  this  is  done  in  Ref  8),,  then 
the  resulting  form  of  the  expressions  for  61  and  6G  is  such  a3  to  provide  an 
easier  route  to  our  6b  '(with  <  ,  rather  that  k  ,,  effectively  involved  from  the 
outset)'.  For  both  6r  and  6w  ,,  however,  the  approach  via  Lagrange's  form  of 
the  equations  is  much  simpler. 
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6.3  The  perturbation  5w 

The  analysis  for  6w  is  much  more  like  the  5r  analysis  than  the  5b 
analysis,  because  each  u£  can  again  be  treated  separately  throughout.  There 
are  two  complications,,  however.  First,  ,(122)  effectively  involves  cos  2v  and 
sin  2v  f/  not  just  cos  v  and  sin  v  (we  see  this  at  equation  (143),,  following), 
and  this  means  that  the  values  J  -  -k  ±  2  are  special  as  well  as  j  -  -k  t  1  . 
Second,  we  cannot  take  5w  to  be  zero  for  any  of  these  special  cases,,  since  the 
constants  in  6e  and  6M  must  now  be  assumed  to  have  been  already  assigned, 
formula  for  the  four  special  5w  will  be  obtained  in  section  7.4.  Actually,,  a 
fifth  special  case  emerges,  corresponding  to  J  -  -k  and  a  zero  denominator 
k  ♦  J  ;  Sw  for  this  case  can  be  set  to  zero,  since  we  still  have  (for  each  k) 
the  constant  in  5u  ,,  as  yet  unassigned,,  available  for  the  purpose  -  the 
constants  for  5u  are  determined  in  section  7.5. 

We  start  by  rewriting  (122),  as 

5w  *  2q’2  'Ue  sin  v  ♦  eq"1  6M  cos  v). 

♦  ieq"2  ( 5e  sin  2v  ♦  eq-1  5M  cos  2v}  ♦  f  e2q~3  ♦  q-1  ,  (143) 

where  6e  ,,  6M  and  51  are  available  from  the  integrals  of  (76)',,  (105)'  and 
(106)'.  The  integrals  for  5e  and  5M  combine  in  a  very  natural  way  and  we  get 


?q'2  (5e  sin  v  ♦  eq"1  5M  cos  v)  -  4q"2  Aik  l  B^.  (e 


Ik  4 

1  -  2k  1  ♦  l  *  2k 


(1  ♦  t  -  k 
k"'+  J  ‘  ^ 


l  k  ♦  J  -  1 


J  ♦  1 


,  f,  1  -  1  -  k  s  1 

3  k  ♦  j  k  ♦  J  ♦  2}  SJ*1J  *- 


'(144) 


sJ-2 


Atk  4  * 

fl*t+k1*l-k 


k  ♦  J  ♦  1 


k  ♦  j  -  1 


sj>,  *  3e 


,k  ♦  J  ♦  2  k  ♦  J  -  2 
1  -  l  -  k 


j  T2J* 


(145) 
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T  *V3  «M  ■  i  eq-2  *lk  I  B4J  {e  I  \  Sj_2 


1  ♦  l  -  2k 


-  R^fVi  +  6e FTT Sj 


q"1  SL  -  (1  -  2«;  A,k  I  B 


■tk  1  kTJ  J 


'msj 


(147) 


We  substitute  the  last  four  results  in  (143)  and  at  the  same  time  (as  in 
the  analysis  for  6r  }  change  the  interpretation  of  J  so  that  we  can  use  the 
same  Sj  in  each  term.  This  leads  to 


-2  a  ?  I*. .of  1  -  fc  ♦  k  1  ♦  t  -  kln 

Alk  i  | 3*  k  *  J  *  2  k  *  J  JBi. 


■k 


*  29 11  k  ~  r  6  k  >~  V  *  3 
1  *  i  -  kl, 


8*2 

1  *  t  -  2k 
k  »  ■ 

1  ♦  4,  ♦  k 


1  ♦  l  ♦  2k 


*  2  k  ♦  J  -  1  p.J+1  *  |®  k  ♦  JTT  *  8  'ITT  j  ♦  i 
.  2  4(1  -  2t)'  *  q2(5  -  t)  ,  g  1  .  I  -  2k  ,  02  1  *  A  -  k 


J  *  1 


2o|2^44^.  3 144424. 6 


k  *  J  «  1 


k  *  J 


9  k  »  J  -  2  J  eiJ 

1  -  t  -  k  1  *  i  -  2k  ] 

k  *  J  -  1  —  -  2  J  “t.J-l 


->fl  ♦  t  ♦  k  1  -  t  -  k  \ 

39  k  *  J  k  *  J  -  2 j0»,J-2 /  SJ  ' 


(148) 


Though  the  algebra  is  tedious,  we  can  now  eliminate  and  ’ 

by  the  appropriate  versions  of  '(54).  If  we  express  the  result  as 

JV  ■  k-2  Atk  z  WJt1  BliJt)  .  VJ>0  By  *  Vj  .,  Byj.,)  Sj  ,  ,(149)' 


the  formulae  for  Vjj  ,,  V^q  and  Vj^  are  initially  very  complicated.  For 
Vjj  ,  in  particular,,  we  start  with 


3Uk  zJ-Lll 

♦  £ '*  i)(k  ♦  j  ♦  2) 


k  -  n  ♦  2  +  6(k  -  l  *  D~ 
k  ♦  J  ♦  2  k  ♦  J  ♦  1 


31(k  -  t  -  D 
(j  ♦  i  ♦  D(k  ♦  J) 


3k  2(k  -  l  -  1)1 

k  ♦  J  k  ♦  j  -  1  J  5 


TR  89022 


45 


is  symmetrically  related  to  this,  but  Vj>0  is  a  great  deal  more 
complicated.  All  three  formulae  can  be  greatly  simplified,  however;  for  Vjj0 
this  was  done  by  a  technique  akin  to  partial  fractions.  The  resulting 
expressions  are 


_ 2 _ § ♦  _L_  .  l 

k  *  J  *  2  k  ♦  J  ♦  1  k  *  j  k  *  j  -  1 


vj,0  ■  8 


l  *  2k  *  1  21  -  1  t  -  2k  ♦  1 

k  ♦  J  ♦  1  "  k  *  J  k  .  J  -  1 


.  ,A %  US.  t  L  .  2U  *  -U  * 
“  k  »  J  ♦  2  k  ♦  J 


tt  *  O  .  i—.-H-i—1- 
k  ♦  J  k  *  J  -  2 


'(151)'  | 


VJ.-1  "  2e(*  “  J)(k  ♦  J  ♦  1  *  k4j  “  k  ♦  J  -  1  +  k  ♦  J  -  2)  '  (152)'  I 

As  a  result  of  this  remarkable  simplification,  it  will  be  observed  that 

Vj  j  j and  vjf.t  Bjfj.j  ,,  in  (149),  have  been  expressed  in  a  very  suitable 
form  for  the  application  of  (56)'  and  ,(60)\  with  l  replaced  by  i  *  1  in  both 

relations,,  to  eliminate  and  respectively,  in  favour  of  By 

(already  present  in  (149))'  and  Bl.1  ^  .  Thus,,  if  we  now  write 

JW  -  +»ll(  l  (W*(0  By  ‘  Wy,  Bj,.,  ,j)Sj  ,  (153) 


Mt.O  ‘  2[k  ♦  J  .  2  '  2  k  *  j  *  k  *  J  -  2| 


sisio:  1 


1  4  6 

k  ♦  J  *  2  '  k  .  J  ♦  1  *  k  *  j 


k  *  J  -  1  k  ♦  J  -  2 


TT  ■  J'55>  I 
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The  final  result  we  require  follows  from  the  substitution  of  (154)  and 
055)  into  (153).  Writing  Sj  in  full,  we  get 

5Mik  •  Aik  |  Tk"']"  2>'fkV  j)Tk  ♦  r-  2>:{C2(1  *  °  ■ k(k  * J)1  BIJ  I 

-  Tk  TJT%TiTT  j  -  o  «"  (XU-  *  jv)'  .  0  56)'  | 

Equation  ( 1 56)  is  the  general  formula  for  dw  due  to  uj[  .  As  with  (130 
and  ( 1 40) •  for  6 r  and  <Sb  respectively,  it  applies  for  all  til;  like 
(140)  but  unlike  (13D»,  on  the  other  hand,  values  of  |j|  up  to  1-1  are 
required  to  cover  all  the  non-zero  terms.  For  each  k  ,  zero  denominators  exist 

for  five  different  values  of  J  :  for  four  of  these  values  (  J  -  -k  ±  1  and 

J  •  -k  ±  2  )',  special  formulae  are  required,  in  place  of  ’(156)t/  as  already  noted; 
only  for  the  fifth  value  ”(  J  -  -k  )'  can  a  term  (for  each  k  )'  be  actually 
excluded. 

Before  proceeding  to  a  pair  of  formulae  for  Niw  K  the  total  number  of 
terras  required  to  express  dw  tor  a  given  value  of  1  ,  we  note  (and  make 
allowance  for)  one  specific  null  term  that  arises  for  each  even  value  of  1  . 

For  k  -  2  and  J  -  1  -  1  ,  we  see  from  (156)  that  the  coefficient  of  B^j  is 
identically  zero  (l.e.  Independently  of  ft  )'.  But  Bj.,  j  is  itself  zero  when 
J  -  ft  -  1  ,,  so  this  specific  term  of  $w^  2  always  vanishes.  Proceeding  to 
N^w  ,  we  first  note  that  the  number  of  J  values  for  each  k  (regardless  of  any 
exclusion)'  is  2ft  -  1  (for  all  ft  i  1  -)'. 

If  ft  is  odd,  then  a  priori  the  value  of  N^w  is  *  D(2ft  -  1)'. 

There  is  one  excluded  value  of  j  for  each  k  *  ft  so  there  are  $-(ft  -  1) 

exclusions  altogether.  It  follows  that 

Nl„  -  *2  ,  ;<157! 

and  values  for  ft  up  to  15  are  given  in  Table  6. 

If  ft  Is  even,  the  a-priori  value  of  Niw  is  -J-U  +  2)(2ft  -  1)  .  There  is 
again  an  excluded  J  for  each  k  *  ft  ,y  amounting  to  H  basic  exclusions ,y  but 
there  are  now  two  other  sources  of  discounted  terms.  We  have  Just  remarked  on 
the  particular  zero  term  that  arises  for  k  -  2  ;  wa  might  prefer  to  allow  zero 
actually  to  be  computed  in  a  general  computer  prograu,  but  here  we  regard  this 
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term  as  an  exclusion.  The  other  source  of  discounted  terms  consists  of  the 
S.  -  1  duplications  that  occur  when  k  -  0  (see  the  footnote  of  section  6.1)'. 
Thus  the  number  of  effective  exclusions  Is  $1  ,  and  from  subtraction  we  obtain 


I  ,  »»„  -  -  1  •  058) 

j 

i  Values  for  l  up  to  16  are  give  in  Table  6. 

|  6.H  Universality  of  results  (non-elliptio  orbits) 

j  Equations  Jl3D  .  ( 1  ’  0)  and  (156)  give,  on  summing  over  k  or  <  as 

5  appropriate,,  general  formulae  for  the  perturbations  fir  ,,  fib  and  fiw  , 

j  respectively,,  due  to  .  It  is  being  tacitly  assumed,  in  the  rest  of  the 

\  Report,  that  we  are  only  considering  elliptic  orbits.  It  is  worth  remarking 

i  here,;  therefore,;  that  (as  follows  by  a  continuity  argument)  the  formulae  are  also 

j  valid  for  parabolic  and  hyperbolic  orbits.  The  formulae  are  effectively 

universal1^,,  in  other  words,  though  they  Inevitably  fail  for  rectilinear  orbits 
(with  infinities  arising  from  zero  p)'. 

7  THE  SPECIAL  CASES,  AND  INTEGRATION  CONSTANTS 

The  main  results  in  this  section,;  obtained  in  section  7.^,  are  the  formulae 
required  to  supplement  (156),,  the  general  formula  for  fiw  .  These  formulae, 
covering  the  cases  J  -  -k  ±  1  and  -k  t  2  ,,  are  forced  by  the  ’constants'  for 
fie  and  5M  ,,  which  are  determined  so  that  certain  terms  (those  for  J  -  -k  ±  1)] 
can  be  excluded  from  fir  .  Though  we  have  omitted  (in  section  *0  the  full 
expressions  for  the  short-period  perturbations,  fic  in  the  elements,,  we  give 
here  the  adopted  'constants'  for  all  the  5  .  Five  of  the  elements  have 
.  constants  chosen  to  suit  fir  ,  fib  and  fiw  ;  for  completeness,  we  start  with 

the  semi -major  axis,,  for  which  the  constants  are  mandated  by  the  use  of  a' 

|  as  3  . 

j  7.1  Mandatory  constants  for  fia 

i  We  go  back  to  the  original  expression  for  fia  due  to  uj[  ,,  viz  (63).  We 

|  can  expand  the  complete  factor  (p/r)**1  in  terms  of  the  , j  the 

expansion  via  the  B?j  in  (65)).  On  taking  Just  the  terra  of  the  expansion  with 
j  -  -k  ,  we  isolate  the  constant  term  that  (for  each  k  ,,  and  a  given  )'  is 
[  mandated  by  taking  "S  -  a'  . 
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The  result  can  be  written  in  the  form  (for  the  ’constant'  component  of 

5alk  $ 


Silk(e)  "  ‘2a<<"2  Atk  sl*2,k  503  k“'  •  (,59> 

7.2  Constants  for  6e  and  6M 

The  task  in  this  section  is  to  derive  the  formulae  for  fieak(0y  and 
5Mfck(c)  that  Initialize  our  taking  the  terms  in  6rlk  for  j  -  -k  ♦  1  and 
-k  -  1  to  be  zero.  These  'constants'  will  complete  the  formulae,,  for  5e  and 
6M  ,,  given  by  the  integrals  of  (76)  and  (105)  respectively. 

We  start  by  observing  that  03D.,  the  general  formula  for  6rik  was 
obtained  by  combining  the  two  different  denominators  from  (129)  and  (130).  If  we 
do  not  combine  the  denominators,,  we  can  rewrite  the  formula  as 

irtk  ■  -ill  -  Dp  Atk  l  (y-pj  .-j  -  ^  j.-.-j]  Bt.u  Cj  ■  0  60)  I 

The  first  denominator  here  is  associated  with  the  I.  summation  of 
section  6.1.  If  this  summation  still  applied  for  J  -  -k  ♦  1  ,  then  the  result 
would  be  an  infinite  coefficient  of  _k+1  C_k+1  .  We  actually  want  this 

coefficient  to  be  ~$-(l  -  1)p  Aik  ,,  since  it  will  then  neutralize  the 
coefficient K  Kl  -  Dp  Aik  ,,  that  arises  without  difficulty  from  the  second 
denominator  in  (160).  The  situation  is  similar  when  J  -  -k  -  1  and  we  want  the 
coefficient  of  C.k.1  ,  from  the  second  term  of  (160)',,  to  be 

- -tU  -  Dp  A^k  (and  not  infinity)]  to  neutralize  the  first  term,  (it  is  recalled 
that  infinite  coefficients  are  avoided,  simply  because  we  deal  separately,;  in 
section  5,  with  the  relevant  terms  of  de/dv  and  dM/dv  .)  What  we  do,, 
therefore,  is  to  obtain  the  coefficients  of  C-k+1  and  C_k-1  that  would  apply 
in  the  absence  of  the  constants  ‘5eg,k(c)  and  5Mi,k(c)'  5  we  can  then  derive  the 
appropriate  values  of  these  constants  to  cancel  these  putative  coefficients. 

So  what  would  the  first-denominator  coefficient  of  Cj  be,,  with 
j  -  -k  ♦  1  ,  in  the  absence  of  the  constants?  There  would  then  be  no 
contribution  from  equation  (123),  but  still  a  cont  ibutlon  from  the  complementary 
(124),  given  by  <$a  .  Its  value  may  be  obtained  from  the  first  term  of  each  pair 
in  (124)  -  the  second  term  does  not  apply  because  it  feeds  separately  into  the 
second-denominator  coefficient  of  C_k>1  which  behaves  normally  as  we  have  seen. 
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But  equation  024)  was  written  down  before  the  ’rearrangement’,  from  (127) 
to  (128) ,,  in  which  the  use  of  J  changed.  This  change  affects  the  B  subscripts,, 
and  it  may  be  seen  that  the  required  first-denominator  coefficient  of  C_k+1  is 

-  a  Alk  (eBlt.k*z  *  *  3eBj,_k),  . 

The  (normally  behaved)  second-denominator  coefficient,,  on  the  other  hand,, 
may  be  written 

£  a  Aik  U  -  1)'  q2  Bji-i,-k*l  . 

To  cancel  the  combined  coefficient  by  use  of  6e^k(cy  and  $Mj,k(c)*- 
let  us  suppose  that 

6etk(c)  "  Aik  x  c*k 
and 

«<UC(0>:  ■  Alk  •*1  Q  »  s-n  <’«2> 

where  x  and  y  are  quantities  to  be  determined.  On  combining  these  for  a 
contribution  to  fir  (cf  (120)),,  we  get  a  coefficient  of  C.k+1  given  by 

-  ±  a  Aik  (x  ♦  y)  ,, 

so  that  one  equation  to  be  satisfied  by  x  and  y 

2(x  *  y)  ♦  eBti.k<2  ♦  ‘tBj,_k,t  *  3eSt,-k  -  .(t  -  1)  q2  ,-k*i  -  0  •  S’63) 

The  complementary  contribution  to  fir  from  (160  and  (162)  leads  to  a  term  in 
C_k_i  .,  of  coefficient 

*  I  a  Aik  ;(x  -  y)  ,, 

and  this  combines  with  two  other  coefficients  of  C_k_1  f/  obtained  as  in  the  last 
paragraph;  the  result  Is  another  equation  in  x  and  y  , 

2(x  -  y)  +  «Bl#.k_2  ♦  ^B4|.k-1  *  3eB£(-k  -  (i  -  O  Q2  , -k~ 1  *  0  •  i164) 
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Solution  of 

063) 

and  '(16-4) 

gives 

x  - 

-  Oq2  (Bt. 

-1,-k+l  *  Bi- 

t ,-k-1 l  ' 

6eBt,-k 

- 

4(Bl,-k»1  ' 

"  Bl,-k-l>  - 

e(Bl,-k*2 

*  B£,-k-2 

and 

y  - 

-  Uq*  :<Bt. 

-1,-k*1  *  b£- 

l.-k-P 

- 

'  Bl,-k-l)-  * 

e(Bi,-k*2 

"  Bfc,-k-2' 

To  get  the  formulae  for 

ieik(c) 

and  6Mtk(o) 

that  we 

require,,  : 

substitute  (165)  and  (166)  into  (161)  and  (162)  .  In  doing  this,,  we  make  two 
siraplif icationss  we  eliminate  B^_1  and  B^_j  by  use  of  (60),  and 

(56),,  respectively  (with  l  replaced  by  l  -  1  in  each  case);  and  we  write 
Bj  k  etc  rather  than  _k  . 

Finally,,  then,,  we  have 

ielK(c),  "  '  *  A£k  '  U  ‘  K  * 

-  (t  -  k  -  *)Bl>k.,  *  eBj,k_2!  00a  ku'  "(167) 


SHtk(c)  ‘  *  6''  <1  *tk  leBt,k*2  -  U  *•  k  - 

-  2keBlk  ♦  It  *  k  -  t0Ba  k_,  -  eBl(k.2]  aln  kM,  _  ( 1 68> 

The  formulae  could,  of  course,,  be  reduced  to  a  smaller  number  of  terms,  by  use  of 
the  fixed-i  recurrence  relation,  '(5*0,  but  the  coefficients  would  then  be  much 
more  awkward;  no  genuinely  simpler  versions  of  067),  and  (168)  have  been  found. 

7.3  Constants  for  ii  and  an 

In  this  section  we  derive  formulae  for  6itk(c)  and  <5nik(c)'  to 
legitim* ze  our  taking  the  terms  for  j  -  ♦  1  and  -  1  in  (IHO),  the 

general  expression  for  $b^K  ,  to  be  zero.  The  analysis  Is  somewhat  simpler  than 
that  in  the  preceding  section,  in  spite  of  the  complexity  entailed  by  •‘he  need  to 
work  with  both  k  and  <  . 
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As  with  6rak  ,,  we  start  by  observing  that  (ISO)  was  obtained  by  combining 
two  denominators,,  which  appear  separately  in  the  preceding  (139).  When 
j  -  -k  ♦  1  ,,  the  second  denominator  becomes  zero  and  no  longer  operates;  from 
the  first  alone  we  get,,  as  the  effective  term  in  (ISO),,  H  At<  Bfc  ^  . 

When  J  »  -tc  -  1  ,  similarly,  the  first  denominator  in  (139)  does  not  operate, 
and  (ISO)  effectively  reduces  to  -U  Ai(C  j  c-k-1  •  These  terras  have  to  be 

cancelled  by  the  use  of  5ij,|<(c)'  and  6Qlk(c)'  ►,  with  appropriate  k  ,,  so  we 
suppose  that 

“Uc<cJ  *  *CH<  t,69)' 

and 

««4k(0)  -  f  S-k  • 

(In  the  last  section  we  were  able  to  include  A^k  in  the  corresponding 
expressions,  (161)'  and  (162),,  but  there  is  no  common  factor  available  now.)'  On 
combining  ( 1 69)  and  (170)  for  a  contribution  to  6b  (cf  '(121))',,  we  get 

±(x  -  y),  C^J1  ♦  ±(x  ♦  y)‘  cbj1. 

If  we  postpone  consideration  of  any  difficulties  associated  with  the 
extreme  values  of  k  that  arise,;  then  we  get  a  pair  of  equations  for  x  and  y  ,, 
on  identifying  *  with  k  ♦  1  and  k  -  1  ,:  respectively,,  in  the  coefficients  of 
and  that  have  been  recorded.  Thus  the  equations  are 

2(x  -  y)'  ♦  l  Al(k+1  B4p.k  -  0  (171) 

and 

2(x  ♦  y»  ♦  l  Atjk-1  Bt|-k  -  0  .  (172) 


Solution  gives 


x  •  ~H(Ai|k*i  *  Ai,,k-1 


(173) 


and 


y 


*l(A, 


'  kt, k-1)Bi,-k 


(17*0 
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The  expressions  in  brackets  can  be  replaced  by  simpler  expressions,  following 
(27)  and  (20),  and  we  also  replace  _k  by  Bik  .  Then  substitution  in  069) 
and  (170)  gives  us  the  formulae  we  require.  Thus 

5ltk(c)  ‘  *Aik  Blk  003  k"'  (175' 

and 

Si!lk(o>  -  ^cs'2  *tk  Bik  3in  k“’  •  (176) 

We  now  have  to  show  that  075)  and  076)  are  valid  for  the  extreme  values 
of  k  as  well  as  'general  values’.  For  k  •  Jl  ,  the  maximum  value  of  k  ,  there 
is  no  difficulty:  071)  involves  a  value  of  l  ♦  1  for  <  in  5b  ,  but  this  is 
all  right  as  BU  -0  (cf  section  6.2);  thus  x  -  y  -  0  for  k  -  l  .  For  the 
minimum  value  of  k  (0  or  1,-  according  to  whether  i.  is  even  or  odd),  it  is  a 
little  more  complicated:  we  consider  the  two  cases  separately. 

If  l  is  even,  with  k  -  0  Cq  and  Cq1  are  the  same  and  we  cannot 

separate  x  ♦  y  from  x  -  y  in  the  combination  of  069)  and  070).  Thus  there 
is  only  one  equation  to  be  satisfied  (instead  of  both  071)!  and  072))',,  and  it 
involves  only  the  identification  of  <  with  k  ♦  1  (not  also  with  k  -  1  )'. 

The  equation  reduces  to 

x  -  ^  Alf1  Bl|0  ,,  0  77), 


which  is  consistent  with  0  75),,  in  view  of  (29).-  The  value  of  y  is 
indeterminate,  but  this  is  appropriate  for  the  coefficient  of  sin  kw1  ,,  which  is 
itself  zero  when  k  -  0  .  If  l  is  odd  (with  k  •  1  ) ,  on  the  oth*r  hand,  the 
validity  of  075)  and  076),  follows  in  essence  from  the  argument  Just  prior  to 
the  establishment  of  equation  (140),,  and  the  details  are  omitted. 

A  last  point  in  this  section  is  noted  as  no  more  than  a  curiosity.  Whereas 
the  rest  of  5itk  and  5fltk  contain  A^k  and  A^k  ,,  respectively,  as  a  factor 
of  every  terra,,  as  Indicated  by  (79)'  and  (82),  these  factors  are  reversed  in  the 
’constant*  terms,,  as  indicated  by  [(175)'  and  (176).  The  point  was  noted,  for 
l  -  3  .  in  Part  1. 
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7.1*  Forced  terms  in  6w 

We  now  have,  for  each  ,  only  one  ’constant*  at  our  disposal;  denoted  by 
6utk(c)  *  wil1  be  deter!nlne(1  In  section  7.5  so  as  to  validate  the  nulling  of 
the  term  for  J  -  -k  in  the  formula,  (156)’,  for  6wai<  .  For  J  -  -k  ±  1  and 
-k  t  2  ,  on  the  other  hand,  we  are  forced  to  accept  non-null  terms  that  arise, 
via  (122),  from  (167)  and  068),  the  formulae  for  <5eik(c)  and  5Mlk(c)  .  We 
derive  the  formulae  for  these  terms  in  the  present  section,  which  must  therefore 
be  regarded  (from  the  reference  viewpoint)  as  a  completion  of  section  6. 

For  each  of  the  four  speoial  values  of  J  t<  in  principle  we  embark  on  a 
procedure  that  is  similar  to  that  employed  in  section  7.2.  The  basis  of  this 
procedure  is  that  we  re-determine  the  expression  for  fiw  with  the  appropriate 
terms  in  6e  and  <5M  (which  in  their  general  form  would  lead  to  infinities)' 
replaced  by  ’constants*  proportional  to  x  and  y  (as  defined  by  (I6t)'  and 
(162));  the  only  difference  (in  principle)  from  section  7.2  is  that  we  do  not 
have  unknowns  to  solve  for,,  so  that  the  procedure  is  ’direct’.  In  practice, 
however,  because  the  basic  formula  for  5w  ,,  (122),  is  so  much  more  complicated 
than  the  corresponding  formula  for  $r  ,  (120),  it  is  better  to  proceed  a  little 
differently  from  this:  Instead  of  developing  our  four  special  formulae  more  or 
less  ab  initio,  we  start  four  times  from  the  (final)  general  formula  for  fiw  ,, 
,(156) ,,  and  modify  it  each  time  in  the  appropriate  manner. 

We  start  with  J  •  -k  ♦  1  ;  it  will  be  useful  to  have  a  shorthand  for  the 
denominator  that  becomes  zero,  so  we  define 


J 


d 


k  ♦  J  -  1 


(178) 


for  general  J  and  will  eventually  set  d  -  0  .  In  terms  of  d  ,,  we  can  rewrite 
,05H)  and  (155)’,,  which  on  substitution  into  (153)  give  the  final  (156),,  as 


W 


1,0 


i  *  k  ♦  1  i  ♦  1  ^  l  -  k  ♦  1 ' 
d  ♦  3  ,<:ao  d-1 


(179)  | 


wl.-1 


■  2(1  -  1) 
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There  is  no  difficulty  with  (179).  but  in  (180)  a  zero  denominator  appears  when 
d  is  set  to  zero.  We  will  find  that  this  denominator  disappears  when  the 
'replacement  procedure’  is  complete. 

The  terms  we  have  to  replace  come  from  the  first  term  of  (143),,  for  which 
the  general  expression  is  given  by  ’(144).  in  the  latter  equation  the  'offending 
terms'  consist  of  the  first  of  each  pair;  on  changing  the  way  J  is  used,  as 
usual,  the  combination  of  the  three  terms  in  question  is 


Alk 


k  +  J  -  1 


LLLijk  o 
k  ♦  j  -  1 


3e 


1  -  i  -  k 
k  +  j  -  T 


j-'} 


where  the  zero  denominators  are  evident  a3  soon  as  we  set  J  -  -k  ♦  1  .  The 
replacement  terra,  also  based  on  the  f  -st  terra  of  (143),  emerges  when  6e  and 
6M  are  set,,  following  ( 1 61 )‘  and  (162),,  to  Aik  x  C_j<  and  Alk  e^q  y  S.k 
respectively.  The  resulting  term  in  S_k>1  is  A^k  q“2  (x  ♦  y),  S.k*j  (in  the 
analysis  for  J  -  -k  -  1  it  will  be  the  term  in  S-k-1  that  we  need),,  so  the 
required  change  in  fiw  is 

q-2  Atk  fx  *  y  .  +d-1C«<k  -  1  -  1)  Bt>_k»2 

*  2(2k  -1-1)  Bti.k„  ♦  3e(k  *1-1)  Bti.k3l  S.k»,  . 

But  x  ♦  y  is  given  by  (163).,  so  this  can  be  expressed  as 

Aik  |(t  -  1)  q2  Bt.ti.k.,  -  d-'Ce(J  *  t)  Btj_k.2 

♦  *  A  -  1?  Bt|.k„  ♦  3e<J  -  i)  Bt  .k])  S.k„  . 

We  set  J  -  -k  ♦  1  In  this  expression  and  then  invoke  versions  of  (56)  and  (58); 
as  a  result  it  simplifies  to 


±(t  -  1)  Atk  (1  -  2d'1);  Bt.1>.k+,1  S.kM 


Since  only  Bi_1>,k+1  ,,  i.e.  ,,  is  present  in  this  expression,  the 

change  in  6w  can  be  represented  as  >.  change  in  W^  (of  (153)  and  (155))  of 
amount  4(1  -  1)(1  -  2d'1)  ;  thus  ( 1 80)  is  to  be  replaced  by 


W 


1,-1 


(181) 
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It  Is  now  legitimate  to  set  d  -  0  In  (179)  and  (l8l)\  the  result  being 

W*,o  -  ”■§■  (2t  -  k  ♦  2) 
and 

wtf-1  ■  -4U  -  D'  . 

On  substituting  in  (153).  we  get  the  first  of  our  four  special  formulae;  it  can 
be  written  (with  a  change  of  sign  in  the  second  suffix  of  each  B  )' 

iwi.k,-k*1  "  ‘  ^  Ajkl(2t  k  +  2) 

♦  U  -  1),  sin  (kw*  +  v);  .  (182) 

The  second  special  case  is  with  j  -  -k  -  1  .  This  is  the  twin  of  the 
first  case,  so  it  involves 

d  -  k  ♦  j  ♦  1  '(183) 

and  the  expression  for  x  -  y  given  by  (164)'.  There  is  no  point  in  a  virtual 
repetition  of  the  analysis  in  detail,  so  we  proceed  direct  to  our  second  special 
formula;  it  can  be  written 

5wlk,-k-1  "  ^  Aik^2i  *  k  *  2)  0i,,k*1 

♦  U  -  D  sin  (kw*  -  v)  .  (184)' 

The  symmetry  between  '(18?)  and  (184)  is  obvious.  For  k  -  0  ,,  of  course,  the 
equations  reduce  to  the  same  formula,  and  the  footnote  of  section  6.1  is  again 
relevant. 

For  our  third  special  case,  we  require  to  set  J  »  -k  ♦  2  ,  so  we  start  by 
defining 


d  -  k  ♦  J  -  2  . 


Then  (154)  and  ,(155)  take  the  form 


Wi,0 


2  *  *  *  ♦  1  .  2  1  »  1  ♦  1  ~  k  +  1-  j 

K  d  ♦  4  d  ♦  2  d  j 


(185) 


(186) 
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-2(1  -  1), 


4 

d  ♦  3 


4 

d  +  1 


,087). 


with  the  potentially  zero  denominator  in  both  expressions. 

The  terms  of  the  general  (156)  that  we  have  to  replace  now  come  from  the 
second  term  of  0**3)»,  for  which  the  general  expression  is  given  by  (145).  The 
’offending  terms’  are  the  last  three  of  the  six  that  occur  in  the  latter 
equation;  with  the  usual  re-interpretation  of  J  ,  these  three  terms  in 
combination  give 


+  eq-2 


Atk 


i  e 


l  ♦  l  -  k 
k  ♦  J  -  2 


♦  2 


1  *  1  -  2k 
k  *  J  -  2 


*  3®  tttttt  • 


The  replacement  term  is  now  -fr  eq’2  (x  *  y),  *  80  the  required  change  in 

<$w  is 


i  eq*2  Alk  |x  *  y  *  4<J",[e(k  -  t  -  D8t._k.2 

♦  2(2k  -  A  -  1),  8t>.k*,  *  3«<k  *1-1)'  Bj..kj|  S.k,2  . 

With  x  ♦  y  given  by  063).,  this  can  be  expressed  as 


+  eq-2  Atk  (U  -  I)q2  Bt_,  ,_k*,  -  <r'[o(J  *  1  -  1>Bti.k,2 

*  2<2J  ♦  1  -  3)Bt|.k„  *  3e(J  -  1  -  1)8t>.k3|  S.k»2  . 

We  set  J  -  -k  ♦  2  in  this  expression  and  then  invoke  versions  of  (54)', 
(57)  and  (58)  to  eliminate  and  8^  _k  in  favour  of 

Bt-1  -k+2’  Thi8  simplifies  the  result  (for  the  change)'  to 

-  *  Atk(2d-'  -  1)1(1  -  k  *  1  *  d)Bt._k,2  -  U  -  1>81_,i_k.2l  S.k+2  . 

It  now  follows,  from  '(186)  and  '(187),  that  the  altered  values  of  W^  Q  and 
Wj  .j  ,  to  be  substituted  in  (153).,  are  given  by 

wi.o  •  -zHr),u-k  -  1  * d)  <,88> 
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wt,-i  ■  -(1  -  '’(ttt'tti*  drfr'rrr 

v»9  can  now  set  d  •  0  K  getting 


(189) 


W1>0  -  -  ±U  +  k  ♦  5)'  (190) 

and 

wl,-1  "  *  'l  ‘  (191) 

Thus  the  substitution  gives,  for  our  third  formula, 

iutk,-k*2  ‘  *A  Alk  i3<*  *  k  *  5>Blik.j  -  I9U  -  DBt., >k_2}  sin  (ko'  ♦  2v)  . 

...092) 

Our  final  special  case,,  with  J  -  -k  -  2  ,:  Is  the  twin  of  the  preceding 
(third)  case  and  Involves 


d  •  k  ♦  j  *  2  . 


(193) 


We  will  not  go  through  the  analysis  In  detail,.  In  view  of  the  symmetry,,  but 
proceed  directly  to  the  final  formula;  It  Is 


s"tk,-k-2  "  A  Alk  |3U  -  k  ♦  5)8tik<2  -  19(1  -  '>Bt-t,k.jl  sin  (kw'  -  2v)  . 

...(199) 


7.5  Constanta  for 

To  complete  section  7,  it  remains  to  determine  the  constant,  5utk(c )'  *• 
that  legitimizes  our  taking  the  tern  for  j  -  -k  in  (156),  the  general 
expression  for  6w^k  ,  to  be  zero.  We  already  have  6M^(C)  ,  given  by  ,0  68),  so 
we  only  need  to  determine  6lik(c)  ,  the  constant  In  6L^k  ,  for  6u4k(0)'  to  be 
known  at  once. 
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We  start  by  now  defining 


d 


k  ♦  j 


:(195) 


so  that  (15**),  and  (155),  in  their  general  fora  with  potentially  zero 
denominators,  are  rewritten  as 


and 


"*,0 


2 


'l  ♦  k  ♦  1 
d  +  2 


2 


d 


l  -  k  ♦  1  ' 
d  -  2 


-2U  -  ’’(jTT 


4 

d  ♦  1 


6 _  4  _j _ ; 

d  ~  d  -  1  d  -  2 


(196)  | 


(197)' 


The  zero  denominators  will  disappear  when,  for  the  coefficient  of  S.^  in  (153),. 
we  replace  a  quantity  occurring  in  the  general  analysis  with  a  quantity  based  on 
5Mtk(c)  ^eik(c)  not  b0in8  involved);  then  <L^(0y  is  defined  to  null  this 
resulting  coefficient. 

The  quantity  to  be  replaced  derives  from  the  combination  of  (146)  and 
(147).  With  our  usual  re-interpretation  of  J  ,  we  can  write  the  resulting 
coefficient  of  a3  the  sum  of 

t  «<r2  Atk  d*1  |e(t  *  t  -  K)Bl(j.2  *  2(1  *  l  -  2k)BtJ., 

♦6e(l  -  l)Btj  ♦  20  ♦  t  ♦  2k)et>j_,  ♦  o(l  *  l  ♦  k>Bl  J_2} 

and 

(1  -  2t)  Atk  d->  Bjj  , 

associated  with  e2  q"3  $M  and  qw*  5L  respectively.  The  quantity  that  has 
to  replace  this  coefficient  is  available  immediately  from  ( 1 68) ,  but  it  is  more 
convenient  to  back-track  a  little  and  take  it  instead  as  f  eq"2  Aik  y  ,  with  y 
given  by  (166).  Thus  the  replacement  coefficient  may  be  written,,  with  j  rather 
than  -k  in  the  8  subscripts,  as 


-  T  eq‘2  Alk  |e  Bj,j,2  *  4  -  4  Bj,j.,  -  e  B4>j_2 

-  U  -  1)d2  (B(.,  jt)  -  8,-,^.])!  • 
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On  subtracting  from  this  the  coefficient  being  replaced,  we  get 

-  T  Ajk  q*2  <r'  |3e2U  ♦  J  ♦  ♦  6e(t  ♦  2j  * 

-  ZW21  -  1}  ♦  e2( t  -  5)3B4j  ♦  6eU  -  2J  ♦  OBj  j., 

♦  3e2U  -  J  *  DBjj.j  -  3eq2  d(l  -  -  Bj.,^.,)}  . 

By  application  of  versions  of  the  five  relations  (5*0  to  (58) ,,  we  can  eliminate 
Bt,J+2  and  *•  then  and  *  and  finally  and 

Bi,-1,J-1  *  *hi3  reduces  the  foregoing  expression  to 

*  Alk  d-'|C2(t  ♦  1)'  -  3Jd]  By  ♦  6(1  -  DBj.j  j)  , 

which  represents  (when  j  -  -k  )'  the  adjustment  required  to  the  coefficient  of 
S-k  in  6w  .  On  adding  the  appropriate  contributions  to  ( 1 96)'  and  (197)  we  get 

"l.O  ■  Bp'TVT1'  ♦  11  ;  *  i  '■  -  33 )  :(198)  ] 

and 

wt,-i  ■  '2U  •  °(d4r' r^T'r^r*  r^r)  •  ('99)  I 

We  can  now  set  d  -  0  (i.e.  J  -  -k  )\  getting 

Wt>0  -  8k  (200) 

and 

-  0  .  (201) 

These  results  mean  that,  in  the  absence  of  6Lik(c)  •  we  would  have 
kAlkBlk  33  the  coefficient  of  S-j<  in  dw^  ;  so,  to  null  this,  we  take 

^lk(c)'  “  " kcJ  A£k  Bik  3in  kw'  •  (20?) 
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Frco  0  68)  and  the  definition  of  »  ,  we  now  get 

s*tk(o):  *  "  i  e'1  %  ie  Bt,k>2  "  <*•  *  k  '  k)Bljktl 

♦  2keBlk  ♦  U  -  k  -  <DBlk_,  -  e  Bi  k.2)  sin  ku'  .  ,(203) 

Finally,  60lk(o)  is  given  by  0  76),  so  our  desired  formula  for  Soigk(0y  is 

s“tk(c)'  "  "  *  e_1  Alk  ie  Bl,k*2  -  <*•  ♦  *  -  k)B1>k<1 

♦  2kes-2  Blk  *  U  -  k  -  *)Bt  ^  -  e  Btk_2l  sin  km1  .  <20U) 

It  is  remarked  that  es  is  free  of  singularity  (as  would  be 

expected), y  since  Aik  contains  sk  as  a  factor,,  so  that  ks”1  A^k  is  non- 
singular.  We  also  have  the  non-singularity  of  s  $Qj|<(c)'  ,,  given  by  '(176),,  for 
the  same  reason. 

8  RESULTS  EXEMPLIFIED  FOR  t  FROM  0  THROUGH  4 

To  illustrate  the  main  results  of  this  Report,  derived  for  general  l  '(>  0) 
we  use  them  to  derive  results  for  the  particular  cases  l  •  1 1<  2,  3  and  4  .  We 
start  with  an  analysis  for  l  -  0  ,,  a  case  not  covered  by  the  general  formulae  - 
their  failure  for  t  &  0  steins  from  the  fact  that  the  expansion  (33)  is  then 
inherently  infinite,,  and  not  Just  ’effectively’  so  ,(cf  Table  4)'.  Both  the  cases 
%  -  0  and  l  -  1  (analysed  next)  are  actually  trivial,  since  the  ’perturbed 
motion’  can  '(in  each  case)  be  looked  at  from  a  viewpoint  which  makes  it  pure 
Kepler ian  (unperturbed) .  The  interest  in  these  cases  then  lies  in  the 
interpretation  of  the  perturbation  formulae,  which  relate  to  the  nominal  mean 
elements  c  »  in  terms  of  the  ’true’  (fixed)  elements  of  the  effective  Keplerian 
orbit  -  the  elements  of  the  latter  will  be  denoted  by  c-j*  . 

For  i.  -  2  and  3  we  write  down,,  from  the  general  results,,  the  specific 
formulae  for  dr  ,  6b  and  6w  that  were  given  before  in  Refs  1  '(’Part  1’)  and 
2.  Both  these  papers  gave  also  the  specific  c  that  complement  6r  ,,  6b  and 
6w  ,  and  Part  1  gave  the  6c  that  underlie  them  -  the  6c  for  l  •  2  are  well 
known,  having  been  given  by  many  authors.  For  t  -  4  we  summarize  a  complete 
.(first-order)  solution,  giving  the  c  as  well  as  6r  ,  6b  and  6w  ,  the 
coordinate  perturbations  (6r,,  6b,  6w),  like  the  general  formulae  from  which  they 
are  derived,  have  not  been  published  before. 
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8.1  The  triTial  (but  exceptional)  case  t  -  0 


From  (4),  we  have 


U0  -  -  yJ 0/r  .  (205) 

This  is  confirmed  by  (15)*,  in  which  k  is  restricted  to  zero  so  that  Uq  «  ujj  ; 
a0  q  and  Aq(1)  are  both  unity,  so  Aq  q  -  JQ  by  '(14).  Thus  the  effect  of 
J0  is  to  reduce  the  power  of  the  central  force  as  indicated  in  Ref  4,  the  value 
of  the  overall  ’true'  power  being  given  by 

uT  -  u(1  -  J0)  •  (206) 

The  orbit  can  be  fully  represented  by  yT  and  the  'true  elements'  » 
but  it  Is  instructive  to  exhibit  the  behaviour  of  the  osculating  elements  (as 
well  as  the  perturbations  6r  ,  6b  and  6w  )'  relative  to  the  y  originally 
assumed.  As  the  general  results  of  the  paper  do  not  apply  when  i  -  0  ,  it  is 
simplest  to  derive  formulae  from  the  original  planetary  equations  directly. 

There  are  no  out-of-plane  effects,,  even  as  a  'trivial'  phenomenon,  since  we  at 
once  get 


6i  *  6fl  •  6b  •  0 


(207), 


so  that 


T  -  iT  and  -  nT  .  (208) 

From  (62)  and  (205)  it  follows  that 

6a  -  a  -  a'  -  -2JQa2/r  .  (209) 

This  is  a  first-order  relation,,  as  usual,  with  a  on  the  right-hand  side 
interpreted  as  a'  (-  5  );  it  becomes  exact  if  a2  is  Interpreted  as  aa'  . 
There  should  be  no  surprise  that  (osculating)  a  variei  around  the  orbit  (unless 
it  is  circular):  this  results  from  the  use  of  the  'wrong'  y  ;  with  the  'right' 
y  (uT),,  the  osculating  a  would  have  the  fixed  value  aT  .  Since  (with  r^  «  r  ) 


62 


2£0 
r  ' 


(210) 


where  all  three  expressions  identify  with  V2,,  V  being  the  orbital  speed,  It 
follows  that 


Sj  •  a’ (1  “  Jq)  ;  (211) 

this  is  an  exact  relation. 

The  planetary  equation  for  p  t/  [(71).  gives 

P  -  0  (212)' 

and  it  is  found  best  to  take  the  integral  of  this  to  be 

$p  -  p  -  p  -  -2J0a  .  (213) 

Then  (209)  and  (213)  lead  to  the  non-singular  perturbation  for  e  given  by 

Se  -  e  -  5  -  -J0  cos  v  ;  (21  U)' 

it  turns  out  that  we  cannot  get  a  simpler  expression  for  our  eventual  5r  by 
altering  the  implicit  constant  in  (21 H)',  based  on  the  explicit  constant  in  (213)'. 
We  introduce  pT  t/  and  hence  eT  ,  by  noting  that 

UP  -  UTPT  .  (215)' 

since  both  quantities  identify  wltn  h2  ,  where  h  is  the  angular  momentum.  It 
follows  from  this,  using  (206)  and  (213).  that 

PT  •  F  -  J0a(1  •  e2),  .  (216) 


Thence,  using  (211),  we  have 


Cj  •  U{ 1  ♦  Jq)  • 


(217) 
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The  planetary  equation  for  $  ,,  (83),.  leads  to 
■ji  -  -J0  e"1  cos  v  ,, 

and  we  take  the  integral  of  this  to  be 

$ij)  -  ip  -  ip  -  -Jq  e"1  sin  v  K 

since  (it  turns  out)'  we  cannot  get  a  simpler  expression  for  6w  by 
(implicit)  constant.  Hence  also 

■  w  -  w  ■  _Jq  e-1  sin  v  . 

We  introduce  ^  via  Vj  ,  noting  that  (6),  taken  with  and  without 
(and  with  r-  •  r  ),  yields 

(p  -  pT)/r  •  (e  -  eT)'  cos  v  -  e(v  -  vT}  sin  v  . 

This  leads  to 


v  -  yT  -  J0  a-1  sin  v 


and  hence  (taking  u  -  uT) 


u  -  Uf  -  -Jq  e”1  sin  v  . 


From  (220),  we  now  see  that 


Wj  •  w  . 


The  planetary  equation  for  p 


,(93),  leads  to 


ifi. 

dv 


-2J0  q  r/p  ,, 


(218), 

(219) ' 
changing  the 

(220) ' 

T-suffixes 

(221)' 

(222) 

(223) 

(224) 

(225) 
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and  this  introduces  a  difficulty  in  the  analysis  for  M  ,  since  (for  the  first 
time)  we  effectively  have  a  negative  power  of  1  ♦  e  cos  v  .  The  simplest  way  to 
deal  with  the  situation  is  to  change  the  integration  variable*  from  true  anomaly ,, 
v  ,  to  eccentric  anomaly,,  E  .  Instead  of  (225),  we  have 

■  -2Jo  ■  :<22«> 

the  Integral  of  this  is  evidently  secular,,  rather  than  short-periodic,  but  for 
convenience  we  use  the  notation  appropriate  to  short-period  perturbations  and 
write  (with  the  most  useful  integration  constant) 

6 p  -  -2J0E  .  (227)' 

■  Tne  secular  perturbation  in  P  ,,  that  has  Just  emerged,  is  dealt  with  in 
the  usual  way  by  choice  of  a  suitable  value  for  7!  ,,  not  compelled  to  be  equal  to 
n’  .  With 


iVj.2  ^3  .  WT  •  u{1  -  J0)'  ,,  (228) 

we  naturally  take 

n  •  ny  -  n* / ( 1  -  J0)'  (229) 

exactly,,  the  formula  being  compatible  with  (211).  (Equation  (229)  is  in  the 
spirit  of  '(115)*  though  not  Just  a  particular  case  of  the  earlier  equation,  which 
is  only  valid  for  t  >  0  .)]  Then  (209)  and  (229)  give 

•  n  -  n  ♦  J0n(3a/r  -  1>  .  (230) 

and  hence 

/  -  Tit  ♦  Jq  '( 3E  -  M)  .  (23D 


*  The  use  of  E  as  integration  variable  lead3  to  a  (finite)  solution  of  the 
general  problem  when  i  <  0  .  The  analysis  is  more  complicated  now,  however,, 
as  v  has  to  be  replaced  by  E  in  occurrences  of  sin  u  ,  Induced  by  the 
factor  (sin  3)  of  ,  as  well  a3  in  the  basic  (negative)  power  of 

1  ♦  e  cos  v  that  arises. 
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In  combination  with  (227)  we  nay  now  write 

5L  -  J0(E-M)  -  JQ  e  sin  E  (232) 

(using  Kepler’s  equation)'.  Using  (219),  we  finally  get 

$M  -  H  -  H  -  jQe'Vq  sin  v  +  e2  sin  E)  .  (233) 

We  introduce  via  ET  ,  since  from  the  equation 

r  -  a(1  -  e  oos  E)  (23^) 

we  get  (with  r  -  rT  ,,  and  both  a  -  aT  and  e  -  eT  known) 

E  -  Et  -  jQe^q’1  sin  v  .  (235)' 

Then  subtraction  of  the  versions  of  Kepler’s  equation  for  M  and  MT  ,y  with  the 
aid  of  e  -  eT  again,  leads  to 

M  -  Mj  •  JqC"1  fa  8ln  v  +  ®2  sln  E>  *  (236) 


so  that  (233)  yields 


MT  -  H  .  (237)' 

Thus  four  of  the  Cj-  are  the  same  as  the  corresponding  X  •  the  only 
differences  being  for  the  elements  a  and  e  .  Further,  we  can  apply  (120)  and 


(122)  to  the  $t  i 

,,  getting 

$r  »  jQa(e2q"*  sin  v  sin  E  -  1) 

(238) 

and 

fiw  -  J0eq*2  sin  v  [(2  ♦  e  cos  v)  . 

1(239) 

In  view  of  (207)  ,.  this  completes  the  analysis  for  i  -  0  . 
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8.2  The  trivial  ease  1  -  1 

Froo  (4),, 


0,  -  -  uJ,  -Jj-slnS.  ,(240), 

This  transforms  f> 

U,  -  uj  -  -  liA,,,  -ffe  003  U'  (241)' 

by  [(15),'  in  the  general  analysis,  where 


*i,i  -  ji  £ 3  ■.  t2H2) 

by  [(1  4)',,  and  sin  B  -  s  sin  u  -  s  cou  u’  .  As  noted  in  Ref  4,  (240)  implies  that 
£  ♦  Ui  -  p[r2  ♦  z2  -  2rz  cos  (fir  -  B)]'*  ♦  0(J12)  [(243)' 


where  z  •  R  ,  so  the  overall  potential  is  the  same  (to  first  order)  as  for  a 
central  force  towards  the  point  at  distance  z  f/  and  axially  'north' ,,  from  the 
nominal  centre  of  ’unperturbed'  attraction.  (The  precise  representation  of  this 
configuration  requires  that  for  each  t  >  0  J,  has  a  specific  value,,  given  by 
-<-ji  >*• .) 

We  have  three  essentially  equivalent  parometers  (•  ,,  k\  j  and  z  )',  and 

our  formulae  can  be  expressed  in  terms  of  any  one  of  these,  but  it  is  more 
convenient  to  use  a  fourth  parameter,,  X  ,,  defined  by 

X  -  z/p  -  -s”1  Ai  t \  -  -J\  R/p  .  (244) 

Then  the  general  formulae,,  taken  with  t  -  1  ,:  lead  to  the  following:  031) 
gives 


$r  -  0  ; 

(140),,  with  <  -  0  and  J  -  0  ,  gives 

6b  -  Ac 


(245) 


(246) 
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and  (182),,  with  k  •  1  ,  gives 

6w  -  -Xs  cos  u  .  (247), 

Though  (247)  is  taken  from  a  special-case  formula  (because  J  -  -k  ♦  1  )' ,,  the 
general  formula,  .(156),  actually  gives  the  same  result  because  the  second  term  of 
this  formula,  which  is  responsible  for  the  zero  denominator,  does  not  arise.  Our 
formulae  are  consistent  with  (132),,  (142)  and  (157),  which  give 

(Ni,r.  N1>b,  Ni|W)  -  (0,  1,,  O'  ,  (248) 

as  seen  also  from  Table  6. 


Expressions  for  the  6c  can  be  written  down  easily  enough,,  following  the 
general  analysis,  but  it  is  of  more  interest  to  obtain  formulae  relating  the  c 
(mean  elements  relative  to  the  nominal  attraction  centre)'  to  the  c?  (unchanging 
osculating  elements  relative  to  the  ’true’  centre)'.  This  can  be  done  via  the  c 
(varying  osculating  elements)  and  derived  quantities,  since  no  conventional 
definitions  are  involved  in  relating  the  c  to  the  Cf  • 

We  start  with  (243),  which  may  be  taken  to  express  y/rT  •  It  leads  to 

r  -  rT  -  Xp  sin  $  ,,  (249) 

so  that,,  in  view  of  (245),, 

r  -  rT  ♦  Xps  sin  u  .  (250) 

Now  this  is  true  for  all  u  ;  but  3  -  aT  3  -  eT  and  H  -  MT  must  all  be 
independent  of  u  ,  whilst  defining  r  -  rT  via  ,(120).  This  is  only  possible  if 

3  -  aT  ,,  :(250 

3  •  eT  -  Xq2s  sin  «  (252) 

and 

R  »  Mj.  ♦  Xe“1q33  cos  u  .  (253) 


TR  89022 


68 


Thus  we  have  established  three  of  the  desired  relationships;  (251)  could 
also  be  obtained,  more  directly,,  by  identifying  two  different  expressions  for  V2 
(cf  (210)). 

We  can  proceed  in  a  similar  way  to  get  the  relationships  for  i  and  fl  . 
Some  geometrical  visualization  is  needed,  and  we  may  regard  the  difference 
between  b  and  b-j.  as  validly  defined,  independently  of  the  precise  location  of 
the  'mean  orbital  plane'  which  is  involved  in  defining  the  coordinate  b  .  This 
difference  is  given  by  a  projection  of  the  displacement  z  perpendicular  to  the 
(mean)  orbital  plane,,  such  that 

b  -  b-j.  -  cz/r  -  Xcp/r  .  (25*1) 


Then  from  (246)',, 


l>  •  b*  ♦  Xce  cos  v 


;<255) 


(where  f>  is  actually  zero,'  by  definition,  but  this  is  not  relevant  to  the 
argument).  As  with  (250),  this  is  true  for  all  u  ,  whilst  5  -  bT  may  be 
expresed  in  terms  of  T  -  iT  and  H  -  flT  via  (121).  It  follows  that 


and 


T  ■  if  ♦  Xce  sin  o 


(256) 


H  •  Slf  -  Xcs_1e  cos  u  . 


'(257) 


This  only  leaves  the  relationship  between  w  and  to  be  established. 

It  was  not  obvious  how  to  proceed,  analogously  to  (249)  and  (254),,  via  a  formula 
for  w  -  wT  ,  so  the  procedure  adopted  was  based  on  formulae  for  7  -  vT  and 
n  -  uT  .  The  first  of  these  comes  easily  from  ”(252)  and  (253);  thus 

V  -  vT  ♦  4-Xe-1s  {e2  cos  (u  ♦  v)  ♦  4e  cos  u  ♦  (2  ♦  e2)  cos  .*>}  .  (258) 

For  the  other  relationship,  we  need  the  special  formula  '(that  can  be  derived  for 
the  given  geometry) 

c(u  -  Uj.)  •  sin  u  cos  u  (1  -  iT)  -  (1  -  s2  sin2u)(Q  -  flT)  .  (259) 
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From  this,  using  (256)  and  (257)  together  with  the  expressions  (omitted  here)  for 
61  and  6fl  ,  we  get 


u  -  Uf  -  iAs-*  (es2  co3  (u  ♦  v)  ♦  2  cos  u  ♦  e(1  ♦  c2),  cos  w}  .  (260) 

If  we  introduce  also  the  (omitted)  expression  for  5u  ,,  (260)  gives 

n  -  UT  ♦  -J-As"^  (es2  cos  (u  ♦  v)'  ♦  4s2  cos  u  ♦  e(1  ♦  c2)  cos  «}  .  (261)' 

From  (258  and  (261)  we  have,,  finally, 

to  •  wy  •  Ae"1  3”^(s2  *  e2c2),  cos  w  .  (262) 

It  is  worth  remarking,  in  conclusion,  that  (262)  can  be  used  to  infer  the 
formula  for  w  -  wy  that  seemed  less  obvious,  intuitively,,  than  the  formula  for 
b  -  bT  .  We  find  that 

w  -  wj  -  (Ap/rs)'  cos  u  .  (263) 

Now  that  the  missing  formula  is  available,  it  is  much  easier  tc  visualize  its 
geometrical  interpretation,  especially  for  polar  orbits  (s  •  1),  . 

8.3  The  case  l  •  2 

This  time  we  start  by  noting,  from  Table  6  or  the  underlying  formulae,,  that 

(N2,r,  N2,b*  n2,w>  ■  (2*  2*  3)'  '(264) 

so  that  there  are  altogether  seven  ferns  In  the  coordinate  perturbations.  As  in 
previous  papers,,  we  simplify  the  coordinate  expressions  by  using  the  notation  K 
and  h  ,  where 

K  -  l-J2(R/p>2  (265! 

and 

h  -  1  -  |f  .  (266) 
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Then  Tables  1  and  2  give 

a2,0  -  *  to,  A2|1  -  -Kcs  f,  A2,2  “  to  >,  (267), 

and  Table  4  gives 

Bi,o  ■  ®2,0  "  1  k  02,1  "  b2,-1  *  "b  •  (268) 

The  two  terms  of  5r  are  given  by  (130,  with  (k,,  j),  -  (2,  0)  and 

(0,  0)  .  We  get,  immediately,, 

Sr  -  ±Kp(f  cos  2u  -  2h)  ,  (269)' 

confirming  equation  (188)  of  Part  1.  (The  single-term  variable  part  of  this 
formula  has  been  given  by  other  authors,,  of  course}  the  best-known  derivation 
was  probably  that  of  Kozai1?,  but  King-Hele  and  Gilmore  established  the  result 
somewhat  earlier,  in  oquaticn  (A-59)  of  Ref  16.) 

The  two  terms  of  db  are  given  by  {140)'  with  (*,,  J)  -  (1,,  1)'  and 

(1,  -1)'  ,  since  '(1,,  0)1  corresponds  to  an  ’excluded  term’.  We  get 

6b  -  -J-Kecs  (sin  (u  ♦  v)‘  -  3  si n  w}  »,  (270), 

confirming  equation  (189)'  of  Part  1. 

For  6w  ,  we  might  have  expected  five  terms  after  the  exclusion  of 
(k,  J)  -  (0,  0),  .  But  (2,  1)'  is  an  example  of  the  specific  null  term  given  in 
general  by  (2,  t  -  1)]  ,  whilst  the  terms  for  (C,  1)'  and  (0,  -1)  ,  being 
identical,  are  combined.  None  of  the  terms  is  given  by  the  general  formula, 
(156):  the  tern  associated  with  (2,  0)  is  given  by  (192);  the  term  associated 
with  '(2,  -1)  is  given  by  082),,  and  the  pair  of'terms  associated  with  (0,  tl) 
are  each  glvon  by  either*  (182)  or  '0  84).  Overall,  we  get 


6w  •  i^Kjf  sin  2u  ♦  4ef  sin  (u  ♦  «)'  *  3eh  sin  v{  ,,  (271)' 


confirming  equation  (190)]  of  Part  1. 


*  In  writing  down  specific  formulae  for  6r  and  dw  when  i  is  even  and 
k  zero,  we  must  always  remember  to  double  the  coefficient  of  each  term  with 
j  -  0  ,  when  our  intention  is  to  cover  the  corresponding  term  with  -J 
(cf  the  footnote  of  section  6.1). 
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The  long-term  motion,  for  1  «  2  ,  comes  entirely  from  the  secular  rates  of 
change,  ^  and  u  ,  given  (from  section  5)  by  -Knc  and  4Kn(4  -  5f) 
respectively.  They  are  quoted  here,  only  to  remind  the  reader  of  the  additional 
'carry-over*  terms  in  6r  t/  6b  and  6w  that  they  induce.  Ir  Part  1,,  these 
terms  are  included  in  equations  (19*0  -  (196)',  as  opposed  to  equations 
(188)-  (190);  the  intervening  equations,,  (19D,-  093).,  refer  to  the  velocity- 
coordinate  perturbations,,  namely,  5r  66  and  6w  . 

Finally,,  of  course,  since  Jg  for  the  Earth  is  of  order  for  l  >  2  ,, 

the  perturbations  of  order  J22  have  to  be  taken  into  account  for  Earth 
satellites.  Part  1  gives  a  detailed  analysis  of  these  perturbations,,  and  the 
resulting  formulae  constitute  the  principal  results  of  that  Report:  equations 
(320),  (3^3).  and  (359)  of  Ref  1  give  the  contributions  to  5r  ,,  6b  and  6w  , 
respectively,  whilst  the  long-term  effects  are  covered  by  equations  (297>  to 
(309)’. 

8.H  The  case  1-3 
From  Table  6, 


<N3,r.  N3(b. 

H3(W);  -  (5,  6,  9)  , 

'(272) 

so  that  there  are  20  terms,  in  total 

Part  1  we  write 

,  in  the  coordinate  pe-tur bat ions. 

As  in 

H  - 

*J3(IVf>>3  . 

(273) 

Tables  1  and  2  give 

*3,0  -  ±Hc(2  -  5f ) , 

*3,,  •  -  5f)  . 

(27«) 

*3»2  "  *  i Hfc, 

*3,3  “  ^  Hsf  5 

(275) 

also  Table  ^  gives,,  in  addition  to  quantities  we  already  have  from  (268),, 


®3,0  *  1  *  83,1  -  e,  83,2  *  ■  (276) 
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We  start  with  fir  ,  separating  (for  convenience)  the  effects  for  k  -  3 
and  k  -  1  .  For  k  *  3  ,  all  the  a-priori  values  of  J  must  be  included, 

namely,  1,  0  and  -1;  for  k  •  1  ,  on  the  other  hand,  we  exclude  J  -  0  .  Then 

,030  gives,  corresponding  to  the  two  values  of  k  t< 

«r  -  itHpaf  !«e  sin  (3u  *  vj  ♦  15  sin  3u  *  20e  sin  (2u  ♦  u)l  (277) 

and 

fir  -  ffcHpes  (4  -  5f)  {sin  (u  +  v)  -  3  sin  w}  ;  (278), 

these  conform  with  equation  (408)  of  Part  1.  (The  total  fir  is,,  of  course, 
gi\en  by  adding  the  two  contributions.) 

For  fib  ,,  the  effects  are  for  <  -  2  and  k  -  0  ,  and  the  a-priori  values 
of  J  are  the  five  with  |j|  5  2.  For  <  •  2  we  exclude  J  -  -1  ,  and  for 
k  •  0  we  exclude  J  -  ±1  ;  for  <  -  0  we  also  lose  a  term  on  combining*  the 
terms  with  J  «  ±2  .  Then  (140)  gives,  corresponding  to  the  two  values  of  <  ,, 

fib  •  -  ^Hcf  (2e2  cos  2(u  ♦  v)  ♦  15e  cos  (2u  ♦  v) 

♦  20(2  ♦  e2)  cos  2u  -  30e2  cos  2w}  '(279) 

and 

lb  -  -  iiHo(2  -  5f)l«2  oos  2v  -  3(2  *  e2>)  i  ,(280) 

these  conform  with  equation  '(411)'  of  Part  1. 

For  fiw  K  the  effects  are  for  k  -  3  and  k  -  1  ,,  with  the  same  a-priori 

J  values  as  for  fib  .  For  k  -  3  ,  all  five  values  yield  terms,,  but  only  three  of 

them  come  from  the  general  (156);  for  J  -  -1  we  use  (192),  and  for  J  -  -2  we 

use  (182).  For  k  ■  1  ,,  the  term  with  J  •  -1  is  excluded,  and  the  only  general 

te^ra  is  for  J  -  2  ;  the  terms  for  J  -  1  ,  0  and  -2  come  from  (192),  '(1?2) 
and  (184)  respectively.  Corresponding  to  the  two  values  for  k  ,  we  get 


*  In  contradistinction  to  the  previous  footnote,  and  as  noted  in  general  after 
(141)  in  section  6,  the  two  terms  do  not  have  the  same  numerical  coefficient. 
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5w  «  -  -^Haf^e2  co3  (3u  ♦  2v)  ♦  lie  cos  (3u  ♦  v) 

+  4(5  ♦  e2)  cos  3u  ♦  ZSe  cos  (2u  ♦  w)  ♦  50e2  Co*  (u  ♦  2u)}  (281)' 

and 

6w  -  4r -  5f){e2  cos  (u  ♦  2v)  -  2e  cos  (u  ♦  v) 

-  208  ♦  7e2)  cos  u  ♦  9e2  cos  (v  -  u)}  ;  (282) 

these  conform  with  equation  (413)  of  Part  1. 

Expressions  for  the  long-period  rates  of  change  of  the  mean  elements  can  be 
written  down  from  the  formulae  of  section  5.  The  results  agree  with  (from  Ref  1)] 
(373).,  ( 376) *  (38*0.,  (389)  and  (399),  for  3  K  T  ,,  1!  w  and  R  ,,  respectively. 

8.5  The  (new)  case  t  •  4 
From  Table  6, 


Nu.b’-  NM*  "  .01.  11»,  15)’  (283) 

so  that  there  are  altogether  37  terms  in  the  coordinate  perturbations,,  which  we 
obtain  first.  To  simplify  our  expressions,,  we  define 

0  -  lArJn  (s/p)11  .  (28U)' 


Then  Tables  1  and  2  give 

A*f0  -  48 G(8  -  UOf  ♦  35f 2) ,  A*ti  -  480Gcs(4  -  7f)  ,  (285) 

Au(2  “  -320Gf(6  -  7f ) ,,  Ai|>3  -  -1120Gcsf,,  4  -  5600f2  ;  (286) 

also  Table  4  gi/es,  in  addition  to  quantities  we  already  nave  from  (276), 

®4,0  *  1  *  84  ■  ■fe ( 1  +  ie2),  B4t2  ■  lc2»,  ®4,3  “  .(287) 
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We  start  with  6r  ,,  as  usual,  separating  the  effects  for  k  -  4  ,  k  -  2 
and  k  -  0  .  The  a-priori  values  of  j  are  the  five  values  for  which  |j |  52. 
All  values  apply  when  k  »  4  ;  for  k  -  2  we  exclude  J  -  -1  ;  and  for  k  -  0 
we  exclude  j  -  ±1  ,,  whilst  the  terns  for  J  -  ±2  are  identical.  Then  031$ 
gives ,  corresponding  to  the  three  values  of  k  , 

6r  -  -2Gpf2  |6e2  cos  2(2u  ♦  v),  *  35e  cos  '(4u  ♦  v)  +  28(2  ♦  e2)  cos  4u 

♦  105e  cos  (3u  ♦  u)  +  70e2  cos  2(u  ♦  w)}  (288) 

Sr  -  -  8  G  pf  ( 6  -  7f){2e2  cos  2(u  ♦  v)  ♦  15e  cos  (2u  ♦  v$ 

♦  20(2  ♦  e2)  cos  2u  -  30e2  cos  2w}  '(289)' 

and 

6r  -  -  24  G  p(8  -  4 Of  ♦  35f2) Je2  cos  2v  -  3(2  ♦  e2) }  .  (290$ 

(The  dominant  (e-free)  terms  of  (288)  and  (289)  were  originally  given  in  equation 
(A1 08)  of  Ref  16.)' 

For  6b  ,  the  effects  are  for  <  -  3  and  <  -  1  the  a-priori  values  of 
J  being  the  seven  with  |j|  S  3  .  For  <  -  3  we  exclude  J  -  -2  ,  and  for 
k  •  1  we  exclude  J  ■  0  and  J  «  -2  .  Then  [(140)  gives,,  corresponding  to  the 
two  values  of  <  , 

6b  -  -  4  G  csf  {4e3  sin  3(u  ♦  v),  ♦  35e2  sin  [(3u  ♦  2v$ 

♦  28e(4  ♦  e2)  sin  (3u  ♦  v)  ♦  70(2  ♦  3©^)  sin  3u 

♦  1 40e(4  ♦  e2)  sin  (2u  ♦  w)  -  I40e3  sin  3«1  '(291$ 

and 

6b  -  -  4  G  cs ( 4  -  7f)(4e3  sin  (u  ♦  3v$  ♦  45e2  sin  (u  ♦  2v), 

♦  60e(4  ♦  e2)  sin  (u  ♦  v)  -  l80e(4  ♦  e2)  sin  w 

-  20e3  sin  (v  -  «)}  .  (292) 

For  6w  p<  the  effects  are  again  for  k  -  4  t/  2  and  0  ,,  with  the  sane 
a-priori  J  values  as  for  6b  .  For  k  •  4  ,  all  seven  J  values  yield  terns, 
of  which  five  come  from  the  general  (156$,  for  J  -  -2  we  use  (192)  and  for 
J  -  -3  we  use  082).  For  k  •  2  ,  the  term  with  J  -  -2  is  excluded,  whilst 
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for  j  -  3  »,  (156)  gives  another  example  of  a  ’specifically  null’  term  (as  in 
section  8.3);  there  are  non-null  general  terms  for  J  -  2  and  J  -  1  ;  and  the 
terms  for  j  •  0,  -1  and  -3  come  from  092), ,  ( 1 82)  and  (184)  respectively. 
Finally,  for  k  *  0  the  term  with  j  -  0  is  excluded;  the  other  terms  cone  in 
pairs,  being  ’general*  for  J  -  ±3  ,,  from  (192)  and  094)  for  J  •  ±2  ,,  and  from 
(182)  and  (184)'  for  J  -  ±1  .  Corresponding  to  the  three  values  of  k  ,,  we  get 

6w  •  -Gf2|4e3  sin  (4u  ♦  3v)  ♦  31e2  sin  2(2u  ♦  v)  ♦  4e(21  ♦  5e2)  '(*)' 

*  sin  '(4u  ♦  v)  ♦  28(3  4a2)  sin  4u  ♦  28e(7  *  e2)  sin  (3u  ♦  w)' 

+  1 75e2  sin  2(u  ♦  w)  ♦  I40e3  sm  (u  ♦  3«)}  ,  (293)' 

5w  -  4  G  f (6  -  7f)  (2e2  sin  2(u  ♦  v)  ♦  4e(5  ♦  2e2)  sin  (2u  *  v) 

♦  5(8  -  7e2)  sin  2u  -  80e(5  ♦  e2),  sin  (u  ♦  w) 

-  40e3  sin  (v  -  2w) }  |(294) 

and 

6w  -  4  G  e(8  -  4 Of  ♦  35f2){2e2  sin  3v  -  3e  sin  2v 

-  6(24  ♦  5e2)  sin  v}  .  1(295)' 

It  only  remains  to  give  the  expressions  for  the  t  (secular  and  long- 
period)  from  section  5.  They  may  also  be  derived  (as  a  check)'  from  the  author’s 
early  Ref  4;  also,  the  version  of  Kepler's  third  law,,  given  here  as  (303) 
checks  with  equation  05)  of  Ref  3. 

There  are,,  of  course,  no  secular  rates  of  change  in  S  or  T  .  Their 
long-period  rates  are  given  by  ,008)  and  (109),,  with  Just  k  -  2  .  Thus 

t  -  -480  Gneq2f (6  -  7f),  sin  2w  (296) 

and 

f  -  480  Gne2cs(6  -  7f)  sin  2w  .  (297) 

The  secular'  rate  of  change  of  ,,  given  by  (110)  with  k  -  0  ,  is 

TI  -  480  Gnc(4  -  7f)(2  ♦  3e2)  ,  (298) 
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and  the  long-period  rate  (given  just  with  k  -  2  )  is 

t  -  -960  Gne2c(3  -  7f)  cos  2w  .  (299) 

For  the  variation  of  13  ,  we  work  with  i?  ,,  given  by  the  second  term  of 
(111).  Thus  the  secular  rate  is  given  by 

t  -  -120  Gn<8  -  bOf  ♦  35f2)(b  *  3e2)  (300) 

so  that  from  (298)  and  (300)  the  secular  rate  for  13  is 

&  -  -120  Gn{H(l6  -  62f  ♦  49f2)  ♦  9e2(8  -  28f  +  21f2)}  .  (301) 


Similarly,,  the  long-period  rate  for  Tp  is  given  by 

i  -  -2b0  Gnf(6  -  7f)(2  ♦  5a2)  cos  2»  ,  (302)' 

from  which  the  rate  for  U  la  at  once  available. 

finally,,  for  the  variation  of  R  we  deal  with  the  aeoular  perturbation  by 
the  modification  of  Kepler'a  third  law  given  by  "(116).  Thla  glvea 

n2  13  .  u(i  .  288 Gq3(8  -  tor  ♦  35f2)l  ,,  (303) 

based  on  the  perturbation  rate  (residual  to  the  mean  motion) 

ft  -  tbb  Gnq3(8  -  bOf  ♦  35f2)  '(30b) 


given  by  '(113).  Again,  the  long-period  rate,  from  '(113),;  la  given  by 

ft  -  b80  Gnq3f(6  -  7f),  cos  2w  i  (305) 

this  checks  with  (302)  and  the  long-period  rate  for  U  ,,  which  from  (11b)  Is 

t  -  -1680  Gne2qf(6  -  7f)  cos  2u  .  (306) 


lyf.X-- 
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Fop  completeness  tn  implementing  the  (f Irst-order)  effects  of  Jjj  (and  any 
other  Ji  ),,  it  is  necessary  to  incorporate  terms  relating  to  what  Part  1 
describes  as  the  difference  between  mean  and  semi-mean  elements.  The  effects 
induced  by  the  secular  variation  are  included  by  adding  ($/n)n  ,,  (3/n)m  and 
(R/n)m  to  i5  ,  3  and  FT  ,  respectively,  where  $  ,  &  and  H  are  given  by 
(298),  (301)  and  (304),,  and  where  m  -  v  -  M  as  in  Part  1.  The  effects  induced 
by  the  long-period  variation,  on  the  other  hand,  are  allowed  for  via  additional 
terms  in  the  expressions  for  6r  ,,  6b  and  5w  .  Using  (120)  -  (122) ,,  we  find 
that  these  additional  terms  are  given  by 

6r  -  480  Gperaf ( 6  -  7f)  sin  (u  ♦  «)  ,  (307) 

6b  -  240  Ge^ncs  {3(4  -  7f)  cos  (v  -  w)'  -  7f  cos  (u  +  2u)}  (308) 

and 

5w  -  240  Gemf(6  -  7f){e  cos  2u  ♦  4  cos  (u  ♦  «)  -  4e  cos  2u}  .  (309) 

9  CONCLUSIONS 

The  main  function  of  Part  1  of  the  present  trilogy  of  Reports  was  to 
provide  details  of  a  new  theory  of  satellite  motion,  largely  based  on  the  use  of 
a  particular  system  of  spherical-polar  coordinates  in  the  representation  of  the 
short-period  components  of  the  orbital  perturbations.  The  emphasis  was  on  the 
derivation  of  the  second-order  perturbations  due  to  the  zonal  harmonic  Jg  ,  but 
the  first-order  perturbations  due  to  J3  were  derived  as  well.  The  latter 
derivation  has  now  been  extended  to  an  arbitary  zonal  harmonic,  (where  l 
is  positive) ,,  v?ith  the  development  of  general  formulae  of  which  those  for  J3 
were  Just  a  particular  case. 

The  main  formulae,  which  (in  their  generality)]  are  believed  to  be  entirely 
novel,  are  those  for  the  perturbations  in  coordinates.  The  general  terms  of 
these  formulae  are  given  by  the  summations  in  ( 1 31 ) »,  (140)  and  (156),,  for  the 
perturbations  in  r  ,,  b  and  w  ,  respectively.  Terms  that  would  have  a  zero 
denominator  are  excluded  from  these  summations,,  as  a  consequence  of  the  optimal 
definition  of  mean  elements,  except  that  replacement  terms  are  needed  for  the 
perturbations  in  w  ;  the  formulae  for  the  replacement  terms  are  (182),,  ]0  84), 
(192)  and  0  94). 
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The  formulae  for  coordinate  perturbations  are  complemented  by  the  formulae, 
given  in  section  5,  for  the  rates  of  change  of  the  mean  elements..  In  principle,, 
the  integration  of  the  rate-of-change  expressions  is  immediate,  leading  to  the 
secular  and  long-period  perturbations  in  the  elements.  In  practice,,  however,, 
there  are  complications,  as  was  indicated  in  Part  1.  One  of  these  complications 
results  from  the  fact  that  the  expressions  really  arise  as  rates  of  change  with 
respect  to  true  anomaly,  rather  than  time,  and  this  leads  to  additional  effects 
that  are  short-periodic  in  nature.  However,  the  difficulty  can  easily  be  dealt 
with  via  the  concept  of  semi-mean  elements;  the  matter  was  fully  discussed  in 
Part  1,  and  has  been  touched  on  here  in  the  context  of  the  derivation  of  the 
appiopriate  perturbation  terms  for  t  H  (section  8.5).  The  other 
complications  arise  in  the  long-term  evolution  of  the  mean  elements,,  the  chief 
source  of  difficulty  being  the  well-known  singularities  in  the  standard  set  of 
elements.  A  preliminary  consideration  of  these  difficulties  was  included  in 
Part  1,,  but  a  full  analysis  is  held  over  to  Part  3.  which  will  also  give  some 
numerical  results. 

The  main  limitation  of  the  theory  presented  by  the  trilogy  is  apparent  from 
its  overall  title  -  the  gravitational  field  is  assumed  to  be  axi -symmetric,  i.e. 
represented  by  zonal  harmonics  alone.  For  a  complete  field,  with  the  tesaeral 
harmonics  included,  the  author  has  already  published  some  general  formulae  (in 
terms  of  cylindrical  coordinates  rather  than  spherical  coordinates,,  though  that 
is  a  minor  detail),,  but  they  apply  only  to  near-circular  orbits.  The  formulae 
were  originally  given  in  Ref  10,  then  in  Ref  5,  and  finally  as  equations  (92)  - 
;(9U)'  of  Ref  9. 

In  the  formulae  referred  to,  the  inclination  functions  Involve  an 
additional  suffix,  ra  ,,  to  cover  the  longitude-dependent  harmonics.  For  m  -  0  , 
the  functions  reduce  to  the  aJ(1)'  and  Ai(<  of  the  present  paper,  whilst  the 
formulae  themselves  are  then  equivalent  to  truncated  versions  of  the  present 
equations  (130,,  (1*10)  and  (156).  Since  we  now  have  one  set  of  formulae  that 
relate  to  all  inclination  functions,  though  the  formulae  are  truncated  in  regard 
to  eccentricity,  and  another  set  of  formulae  that  are  valid  for  any  eccentricity, 
though  only  relating  to  inclination  functions  for  which  m  •  0  ,,  an  obvious  goal 
is  the  derivation  of  formulae  that  are  ’general'  in  both  respects.  There  is  a 
fundamental  difficulty,,  however,,  arising  from  the  rotation  of  the  gravitating 
primary,,  which  we  are  able  to  neglect  in  the  trilogy  because  it  is  assumed  to 
take  place  about  the  axis  of  symmetry. 
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The  root  of  the  trouble  i3  that  the  disturbing  function  (U)  Is  no  longer 
time- independent  when  the  rotation  of  an  arbitrary  primary  is  allowed  for.  This 
nullifies  our  key  constant,  a'  K  and  leads  ipso  facto  to  the  important 
phenomenon  of  resonance^.  it  will  not  be  easy  to  develop  a  unified  theory  that 
covers  resonant  effects  by  the  same  formulae  as  non-resonant  ones.  4However,  a 
starting  point  is  obviously  the  generation  of  the  formulae  referred  to  (in  the 
preceding  paragraph),  as  being  ’general  in  both  respects';  Appendix  A  gives  an 
outline  of  what  is  involved. 
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Appendix  A 

EXTENSION  TO  THE  GENERAL  GRAVITATIONAL  FIELD 

Extending  the  theory  of  this  Report  to  the  tesseral  harmonics  (sectorial 
Included)  is  an  easier  matter  than  might  have  been  expected,  so  long  as  the 
rotation  of  the  gravitating  body  is  neglected;  i.e.,  we  suppose  the  sidereal 
angle,  v  ,  to  be  fixed.  We  assume  the  potential  to  be  described  by  the  usual 
harmonic  coefficients,  and  SfcB  *,  where  -C^(0  oan  identified  with  the 

zonal  coefficient  Jg,  and  q  is  taken  as  zero.  For  convenience,  we 
introduce  the  polar  equivalents,,  and  ,  where  \  is  longitude  and 

(C ,  Sg£j)  “  (cos  raAg,m,  sin  mA^)  .  (A-l), 

(Note:  Agjjj  is  not  uniquely  defined  if  m  >  1  ,,  and  if  m  -  0  we  set 
Jj,0  -  ,,  so  Jg,#o  raust  be  allowed  to  be  negative.)  It  is  usual,,  in 

practice,,  to  work  with  normalized  versions  of  Cg,n  and  Sg,n  (and  hence  Jgm  )',, 
but  this  is  an  irrelevant  complication  here.  The  potential  due  to  , 
generalizing  equation  '(H)  of  the  main  text,,  is  given  by 

“im  *  £  Jim  P?(sln  6)  cos  «U  -  .  !(A-2) 

The  expansion  of  Ufcn  ,,  in  terms  of  the  orbital  elements,'  is  customarily 
based  on  the  family  of  inclination  functions,  Fgmp(i}  ,  such  that,  generalizing 
equation  (8), 

i 

P™(sin  6)  exp  (iml)  -  l  F£.mp( i )  exp  i{(Jt  -  2p)u  ♦  m(Q  -  v)}  .  ?'-3) 

p-0 

As  already  indicated  in  the  main  text  of  the  Report,,  however,,  we  prefer  to  use 
the  index  k  (•  A  -  2p)  ,,  rather  than  p  .  This  index  only  takes  values  that  are 
of  the  same  parity  as  l  ,,  but  in  the  extension  to  m  >  0  we  have  to  allow 
negative  values  of  k  ,,  so  that  its  range  is  now  from  to  ;  as 
compensation,  we  no  longer  require  the  factor  u^  introduced  at  equation  (8). 
Further,,  we  prefer  the  inclination  functions  to  be  real  for  all  values  of  the 
Indices,  so  we  define,  as  an  unnormalized  equivalent  of  the  ?£m(i)]  in  Ref  9, 

4,(0  .  .“‘k  Ftmpd)  •  (*-*)' 
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The  Ug,m  expansion  will  involve  Q  ,  as  well  as  the  other  elements,  to  reflect 
the  abandonment  of  axial  symmetry,  but  for  convenience  we  work  with  fl*  defined 
by 

fl*  »  Q  -  v  -  ■W  .  (A-5) 

Then  Uj^  can  be  decomposed  into  I  bf  ,  where 

k 

ula  *  f  (B/r)4  Jlm  cos  fku’  *  »<«'  *  Hm'l  i  (A-6); 

this  la  compatible  with  equation  ,(71)  of  Ref  9.  in  which  if  and  X  were  the 
negatives  of  the  present  u'  and  Q'  .  (Compatibility  with  two  other  papers  can 
be  obtained  by  noting  that  F{m  here  is  ii+k  times  the  unnormalized  equivalent 
of  the  Ffcnp  u*,ed  in  Ref  13.  whilst  F£mp  in  Ref  17  is  Identical  with  Fj,np 
introduced  at  (A-3?  here.)' 

Next  /e  introduce  quantities  Aj^k  that  directly  generalize  the  Ag,m  of 
the  max  jxt,  defining 

-  -Jtm  (B/P)»  .  U-n 

We  also  generalize  cj  and  sj  ,,  by  defining 

cjBk  -  cos  Cjv  ♦  ku'  ♦  m(n’  -  x^n)] 

Sj*k  -  'Jin  Cjv  ♦  ku’  ♦  «(fl’  -  X*n)]  ,,  ,(A-8)‘ 

and  henceforth  we  will  omit  the  superfixes. 

Then  (A-6)  to  (A-8)  give 

J  Vr)M  Alnk  C0  ,  :<A-9> 

which  is  a  straight  generalization  of  equation  ’(15)  of  the  main  text.  When 
m  S  k  ,,  we  can  also  generalize  the  preceding  equation  (14)’  by  writing 

Aiok  •  -Jim  WP>1  Hmk  sk'°  '<>  ♦  =  >“  AklnU)  ,,  (A- 10)' 
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where  is  from  the  family  of  ’normalized'  inclination  functions 

Introduced  in  Ref  17  to  generalize  the  Aj(i)  ,,  and  0^  generalizes  aim  ; 
the  formula  for  ctj^  can  be  inferred  from  equation  .(11)  ibid.  It  then  follows 
from  (A-7)  and  (A-10)  that 

Fta'U>  ■  Hmk  3k'n  C  *  e>”  •  JA-")' 

It  was  assumed,  in  (A-10),;  that  m  £  k  .  When  0  £  k  i  m  ,,  a  different 
generalization  became  necessary  in  Ref  17,  leading  to 

‘ink  ■  'Jim  “tak  >m'k  S'  ♦  e)k  aJb(1>'  ,  <A-12)‘ 

where  now  the  formula  for  can  be  Inferred  from  equation  '(13)  ibid}  also 

(A-7)  and  (A-12)  lead  to 

■  “*mk  »"~k  C  ‘  o)k  4«S  •  W-13) 

For  k  •  a  ,  (A-12)  and  (A-13)  are  consistent  with  (A-10)  and  (A-11)>  • 
respectively,  but  otherwise  the  dual  definitions  of  A^m(i)'  and  are 

distinct.  A  further  complication  is  that  an  extension  of  (A-10)  and  (A-11)'  to 
negative  k  is  not  generally  available;  (A-12)  and  (A-13)  still  operate  for 
k  <  0  ,  with  |k|  S  m  ,,  but  there  is  a  marked  lack  of  symmetry  between  the  forms 
of  A£a(i)  and  a£ak  ^or  *<  <  0  in  relation  to  k  >  0  .  The  difficulty  for 
negative  k  is  not  too  serious ,,  however,,  as  F^d )'  can  then  be  derived  from 

Fta<n  ■  (->° 

By  appeal  to  < A- 1 ^ )  as  required,  the  A ink  can  always  be  obtained.  There 
are  advantages  in  the  adoption  of  a  different  ’normalization*  for  the  **««>  ■ 
however,  such  that  they  constitute  a  fully  unified  family  of  functions,  defined 
for  all  k  and  symmetric  in  regard  to  the  sign  of  k  .  The  constants 
must  then  also  be  redefined,  to  preserve  tne  uncnanged,  and  the  connecting 

formu  ae  is  based  on  (A-12)  rather  than  (A-10).  We  intioduce  a  variable  sign 
into  (A-12),  to  make  the  always  positive;  expressed  symmetrically  in 

regard  to  the  sign  of  k  ,  the  connecting  formula  is  then 

*l,»k  -  <->*U‘1<'2)  Jim  (R/P)1  «tak  s“  (7-^]^ 
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The  new  A^U),  are  defined  for  l  i  m  ,  as  before,  but  for  all  k  now, 
not  just  for  |k|  $  i  (quite  apart  from  the  problem  with  k  <  0  )'.  The 
recurrence  relation  for  fixed  m  and  k  is 

(1  -  DU2  -  a2)Aj[jB  -  (21  -  1)|lU  -  Do  -  ak)  *£_,>n 

♦  *.{(*.  -  I)2  -  k2)  aJ.2[|1  -  0  ,  (A-t6) 

which  is  slightly  simpler  than  the  relation  in  Ref  17;  the  starting  values  for 
this  are 


A^U)  -  1  and  A*+,ta(i)'  -  (a  +  Do  -  k  ,  (A-17)' 

though  the  second  of  (A-17)  can  be  dispensed  with  if  we  define  a£_^  m(i)  to  be 
zero.  A  recurrence  relation  for  fixed  l  and  m  is  also  available,,  viz 

(A  -  10(1  *  o)A^’  -  2(n  -  kcH^  »  (A  *  k>(1  -  o)A^’  .  0  ,  ,(a-18) 

in  which  the  symmetry  [(in  regard  to  the  sign  of  k  )  is  obvious;  expressions  for 
Aj[a(i)[  ,,  with  |k|  <  i  can  be  generated  'from  either  end'  by  use  of  Just  one 
starting  value  from  the  pair 

A^U)'  -  CKo  *  1)]*-*  .  (A-19) 

There  is  also  a  recurrence  relation  for  fixed  l  and  k  ,  but  instead  of  giving 
it  we  note  that  the  set  of  15  relations,  each  Involving  A*m(i)'  and  two 
•adjacent*  functions  f^oro  a  three-dimensional  table,  can  all  be  generated  from 
various  subsets  of  Just  three  relations;  one  of  the  simplest  such  subsets 
consists  of  ( A- 1 8 )[  and  the  following  pair  of  relations: 

it  ♦  k)AJf_,  iI0  -  (m  *  AcJA^  ♦  AO  *  c)a£*'  -  0  (A-20) 

and 

2(  A  -  m  ♦  1)A^m.,  *  (m  *  k)(1  -  0)A^o  -  (A  -  k)0  *  c)a£*'  -  0  .  (A-21)' 
(Note  in  proof:  see  Refs  19  and  20  for  computational  aspects  of  these  relations.) 
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We  will  not  enlarge  on  the  advantages  of  the  redefinition  of  A^(i)  and 
aJtmk  »,  but  tw0  disadvantages  must  be  mentioned.  First,,  since  the  factors  In  s 
and  c  in  '(A-15)  can  be  combined  as  ,(!  ♦  c)^n+k^ (1  -  c)4-(m”k)'  ,  we  see  that  a 
negative  power  of  either  1  +  c  or  1  -  c  appears  whenever  m  <  |k[  ,,  and  this 
has  to  be  cancelled  by  a  corresponding  positive  power  that  is  present  in  the  new 
A^U),  .  Second,  the  use  of  (A-16)  and  ( A- 1 7 >  to  compute  A^m ( i >'  under  these 
circumstances  (m  <  |k|),  is  Inefficient,,  since  the  recurrence  process  has  to 
work  through  the  unwanted  functions  with  a  <  l  <  f k |  . 

The  A^m(i)'  are,  like  the  a£(1)'  of  the  main  text,  defined  regardless  of 
parity.  The  constants  apa]<  (and  hence  the  quantities  Ajtnk  )  are  only  defined 
for  l  and  k  of  the  same  parity,,  however,,  and  (as  redefined)  their  only 
property  to  be  stated  here  is  that  of  complete  symmetry,  so  that 

<4.  •  “tm  <  >  °>  ■ 

But,,  Just  as  in  the  main  text,  we  require  another  set  of  constants,  a^nic  ,  and 
quantities,,  A$,n<  ,,  defined  when  l  and  k  are  of  opposite  parity,,  to  allow  the 
formula  for  6b  to  be  expressed.  The  connecting  formula  corresponding  to  (A-15) 
is 

*t»c  ■  M*'**'”  Jta  O'P'*  “Ink  »n  (frl)**  A«(1)  •  I 

The  <n0<  are  available  at  once  from  the  ain<  ,  since  '(cf  (26)  of  the  main 
text,  which,  because  of  the  redefinition,,  is  not  being  directly  generalized) 

•tow  *  «r  •  U'W 

Tables  and  further  properties  of  the  redefined  inclination  functions,,  and  the 
associated  constants,,  will  be  given  in  a  separate  paper. 

By  making  use  of  the  quantities  and  A^m<  we  find  no  difficulty  on 

extending  the  theory,  largely  because  the  treatment  of  (p/r)**1  ,,  in  (A-9).,  via 
the  ,  goes  through  unchanged  from  the  main  text.  Further,,  the  energy-based 
exact  quantity,  a*  ,  is  still  available,,  followirg  the  assumption  that  the 
attracting  body  does  not  rotate.  Thus,,  equations  (65).  (76),,  (88)  ar.d  (105),  for 
6a  ,  de/dv  ,,  ds'/dY  and  dM/dv  ,  -espectlvely,  are  unchanged  apart  from  the 
appearance  of  in  place  of  kfe  .  Equation  (82),  for  dfl/dv  ,  requires  a 
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corresponding  change,  such  that  the  derivative  A^  replaces  A^m  .  This  Just 
leaves  (79),  for  di/dv  ,  for  which  a  slightly  more  explicated  expression  is  now 
required,  to  reflect  the  fact  that  pc2  is  no  longer  an  invariant.  He  have,  in 
fact, 

.  2rapc  Alnk  I  8tJ  Sj  .  ( A- 25) 

Fran  this,  using  the  version  of  dp/dv  corresponding  to  equation  (7*0,,  we  get 

■ji-  -  a-'fko  -  a)  Ajan  l  By  Sj  ;  (A-26) 


in  comparison  with  equation  {79),,  we  see  that  the  only  additional  change  is  the 
replacement  of  kc  by  kc  -  r„  . 

Six  of  the  seven  formulae  that  define  dr  ,  db  and  dw  completely,  for 
...ie  zonal  harmonics,,  are  immediately  applicable  to  the  zonal  harmonics,  so  long 
as  A^k  replaces  A*m  and  the  trigonometric  argument  includes  the  term 
m(Q*  -  Xjtm)  .  These  six  are  031),  and  056),  for  the  general  dr  and  dw  ,,  and 
082),,  084),  092)  and  09*0,  the  four  special  formulae  for  dw  .  In  the 
seventh  formula,  '(140)  for  db  ,  l  A$K  must  be  replaced  by  '(A  ♦  a)hta<  ,  in 
addition  to  the  inclusion  of  the  extra  term  in  the  trigonometric  argument.  (It 
is,,  perhaps.,  suprising  that  the  change  to  (140)  is  as  slight  as  this,  but  it 
would  have  been  even  less  if  Aj,k  and  A^k  had  beer,  defined  to  include  the 
factors  t  and  t  +  m  respectively;  the  reason  for  excluding  these  factors 
was,  essentially,  to  give  a  degree  of  homogeneity  to  (26)  and  (A-24.)' 

Finally,,  of  course,  the  numbers  of  terms  in  dr  ,,  db  and  dw  ,  for  a 
given  J^m  ,,  are  greater  for  m  >  0  than  for  a  -  0  ,  to  reflect  the  distinction 
between  positive  and  negative  k  .  These  numbers  are  otherwise  independent  of 
n  ,  however,  in  consequence  of  which  we  write  the  formulae  as  follows: 


and 


»[r  -  242  -  3i  ♦  1  ,( A-27) 

N^b  -  2A2  -  3t  *  2  (A-18) 


f  2A2  for  odd  A 

|  2 (A2  -  ])  for  even  A 


(A-29) 
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Appendix  B 

THE  QUANTITIES  Btj  ,  AND  B0,i  IN  PARTICULAR 

When  j  2  0  ,,  Bg,j  may  be  expressed  In  terms  of  the  hypergeometric 
function  as: 

BtJ  -  .(*  J  ')  <V2)J  Fp  ~  *  1  .  J— y-  -  j  3  ♦  1  !  e2j  ,,  (B-1)] 

where  (B-1),  applies  for  all  l  .  This  result  is  proved,  In  terms  of  the  function 
Bj[(e)'  in  Appendix  E  of  Ref  4;  it  is  also  quoted,  in  terms  of  the  equivalent 
Hansen  function,  at  equation  (32)  of  Ref  8.  For  (0  S  )  J  <  l  ,,  (B-t),  gives  a 
polynomial  in  e2  ;  for  J  U>  0  it  gives  zero;  and  for  ISO  it  gives  a 
power  series  in  e2  ,,  which  can  be  transformed  into  a  closed  expression  involving 

q  (-  /(I  -  e2))  and  perhaps  8  (-  e/(1  ♦  q))  .  All  this  is  consistent  with 

Table  4. 

Equation  (B-1)  breaks  down  when  J  <  0  ,  owing  to  the  eventual  occurrence 
of  a  zero  denominator  when  the  hypergeometrlc  function  is  expanded.  Since 

BiJ  *  Bi,-J  <B-2> 

however,,  this  break-down  is  of  no  account.  The  Justification  of  (B-2)  comes  from 
the  Hansen-f unction  equivalence  and  the  relation 

■  X-VJ  ..  <B-3) 

which  follows  immediately  from  the  definition  of  Hansen’s  functions  by  equation 
(35)  of  the  main  text. 

The  quantity  Bqj  is  of  particular  interest,  being  the  simplest  of  the 

Bjjj  that  involve  8  .  Once  it  is  known,  the  other  such  Bgj  can  be 

progressively  derived  using  the  recurrence  relations  of  the  main  text. 

From  (B-1)  we  have 


B0, 1  •  -  4«FO,.  t+j  2;  e2>  . 


(B-4) 
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But  as  a  particular  case  of  the  hypergeoaeWc  relation  proved  as  equation  (9.2) 
of  Ref  18,  we  have 

F(t,  H;  2;  a2)  -  2q"H(i,  1i  +;  q2)  -  2F(1 ,,  Hi  Hi  q2)  .  (B-5) 

Also,  it  is  immediate  (from  the  expansion)  that 

F(a,  b;  a;  q2)  -  '(1  -  q2)"b  ,  (B-6) 

independently  of  a  ,,  so  that  (B-5)  gives 

FO,,  14-5  2;  e2)  -  2/q(1  ♦  q>  .  !(B-7) 

Finally, ,  (B-4)  and  (B-7),  give 

®0,1  *  *  Sq’^  *,  (B-8) 

in  conformity  with  the  entry  in  Table  U. 
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Table  1 

THE  FUNCTIONS  a£(1) 


k 

0 

1 

2 

t 

0 

1 

2 

1  -if 

C 

1 

3 

0(1  -  if) 

i  -if 

c 

4 

1  -  5f  *  Vf2 

od  -  in 

i  -  if 

5 

cd  -  7f  ♦  V  r2). 

i  -  if*  Vf2 

od  -  in; 

6 

1  -  V r*  4*f2-Wr3 

o(1  -  if*  Vf2) 

i  -  3f  *  Hr2 

1 


C  1 

i  -  if  c 

c(1  -  V«  1 


-  Hr 


90 

Table  3 

THE  FUNCTIONS  B^<e) 

_ t 


Table  4 


THE  QUANTITIES  8, 


1 

2 

— £eq”7(4  +  e2) 

*eV7 

-4eq*5 

4e2,-5 

-eq“3 

&-3(1  *  2q) 

-Bq'1 

8V 

0 

0 

b 

0 

e 

■j-e(4*  e2) 

W 

9: 


Table  5 

THE  QUANTITIES  E^ 

o 
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Table  6 

THE  NUMBER  OF  TERMS  IN  Sr  ,  Sb  AND  6w 


l 

N*,r 

Ni,b 

1 

0 

1 

1 

2 

2 

2 

3 

3 

5 

6 

9 

4 

11 

11 

15 

5 

18 

19 

25 

6 

28 

28 

35 

7 

39 

40 

49 

8 

53 

53 

63 

9 

68 

69 

81 

10 

86 

86 

99 

11 

105 

106 

121 

12 

127 

127 

143 

13 

150 

151 

169 

14 

176 

176 

195 

15 

203 

204 

225 

16 

233 

233 

255 

I 

I 


9** 


LIST  OF  SYMBOLS 
(Usage  for  the  main  text  only) 

a  semi -major  axis 

?*  energy- based  fixed  value  of  a 

a£(i)  ’normalized*  function  of  inclination 

Afck  quantity,,  based  on  A^(i)  ,  defined  by  (1*0 

Ajk  similar  to  Ajj<  ,,  but  defined  by  ,0?), 

A^  derivative  of  Ajj<  with  respect  to  1 

A^  ks-1  A^  t  c*'  A^ 


b 


BIJ 


a{i 

C 

Cj  (0r  Cj> 
d 

e 


E 


latitude-like  spherical  coordinate  of  (r,,  b,,  w) 
normalized  function  of  eccentricity 
quantity  related  to  B^(e)'  ,  defined  by  (32) 
derivative  of  Btj  with  respect  to  e 
cos  l 

cos  (Jv  ♦  ku’)'  (different  meaning  in  Part  1) 
shorthand  for  k  *  J  -  1  etc  in  sections  7.4  and  7.5 
inclination  function,  quoted  frco  Ref  4 
eccentricity 

eccentric  anomaly  (only  required  in  section  8) 
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rsir’^i}  (R/p)*  (in  section  8,5) 

angular  momentum  (but  1  -  -ff  in  section  8.3) 

4-  J3  (R/p)3  (in  section  8.4) 
inclination 

index  associated  with  multiples  of  v 
zonal  harmonic  coefficient  for  the  Earth 
index  associated  with  multiples  of  u 

|  J2  (R'p)2  (in  section  8.3) 
index  of 

quantity  such  that  L  M  ♦  q<}>  -  n  p 

v  -  M  in  section  8.5  (and  Part  1);  otherwise  an  arbitrary 
integer 

mean  anomaly 


mean  motion 


We  will  require  derivatives  of  the  inclination  functions.  It  is  evident 
frora  (10)  that 


_d_ 

di 


lA^m 


(i  -  k)U  ♦  k  1 ) 
2(k  ♦  1) 


3  A^’.U)  , 


(18) 


from  this  and  X(14)  it  follows  that  the  (partial)  derivative  of  with  respect 

to  i  is  given  by 

■  h  (s/?)*'  “tk  3k'1  {k<=  ❖I).  -  f  Af'U)}  ,(19)  I 

where  f  -  s2  .  Tne  quantity  in  [(curly)'  brackets  is  the  D^(i)'  of  Ref  4.  We 

will  also  require,,  finally,,  the  particular  combinations  of  Aak  and  AjK 
denoted  by  Ajk  and  A^k  F/  and  given  by 

Aik  "  k3-1  Aik  4  C"1  Aik  »  (20) 

the  s*1  and  c”1  factors  do  not  imply  singularities,,  as  they  must  always 
cancel  via  k  Ajik  and  Ajk  respectively. 

The  A  £  ( i )'  and  a^k  (and  hence  the  A^k  ’)'  may  be  computed  with  the  aid  of 
recurrence  relations.  A  fixed  k  was  stipulated  in  Ref  4  for  the  formula 

U  ♦  k)  Aj[(i)  -  (2i  -  1)  c  Ak_1  ( 1 )  -  U  -  k  -  1)  A1J_2(i)  '(21) 

valid  for  l  l  k  ♦  2  with  the  starting  values  Ak(i)  •  1  and  Ak+1(i)  -  c  ; 

(?1)  Is  even  valid  for  l  *  k  +  1  ,  if  an  arbitrary  (but  finite)  Ak-1U)  is 
assumed.  However,-  it  is  usually  more  useful  to  stipulate  a  fixed  l  ,  the 
required  formula  was  giver  by  Merson1 1 ,,  being 


AjCi ) 


c  Afl(i)  - 


U  -  k  -  i)U  ♦  k  ♦  2) 
^(k  ♦  1)(«  ♦  2) 


f  A5^(i) 


(22) 


valid  for  l  -  2  >  k  2  0  with  the  starters  A*(i)  -  1  and  A|"Ui)  -  c  , 

( 22)  is  also  valid  for  k  «  l  -  1  ,,  with  an  arbitrary  (finite)  a|*‘(i)  .  Either 
of  the  two  preceding  'pure'  three-term  recurrence  relations,,  (2i)  or  (22;,  can  be 
usea  with  just  one  'mixed’  such  relation  to  generate  ail  the  relations  connecting 
the  A^(i)  ,  perhaps  tne  simplest  "»ixe>l  relation  (with  neitner  l  nor  k  fixed) 

13 
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LIST  OF  SYMBOLS  (continued) 


Nt,r*  Nt,b*  Nfc,w 
P 

pt<  > 
p*<  >. 

Q 

Qt(e)' 

r 

R 

Rj  (for  R^) 

3 

(or  Sj) 
t 

Tj  (for  Ti^j) 
u 

u’ 


uk 


number’  of  terms  (for  given  t)  in  6r,  6b,  5w 
parameter  (seml-latus  rectum)  of  orbit 
Legendre  polynomial  (of  argument  supplied) 

Legendre  associated  function 
/(I  -  e^) 

normalized  eccentricity  function  quoted  from  Ref  3 
radius -vector  coordinate  of  (r,  b,  w ) 

Earth’s  equatorial  radius 

quantity  defined  by  (125)'  (different  in  Part  1), 
sin  i 

sin  (Jv  +  ku’)  '(different  meaning  in  Part  1)' 
time 

quantity  defined  by  037) 
argument  of  latitude,,  v  +  u 
modifier'  u  (-  v  ♦  «’  »  u  -  i«), 
potential  due  to 
component  of 
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LIST  OF  SYMBOLS  (continued) 

v  true  anomaly 

V  orbital  speed  (used  only  in  section  8) 

V«  Quantities  Introduced  at  (149), 

w  longitude-like  spherical  coordinate  of  (r,,  b,,  w) 

0»  w8,  -1  quantities  introduced  at  (153) 

x,  y  general  unknown  quantities  (different  in  Part  1) 

X*J  generic  Hansen  function  (of  eccentricity), 

oi 

z  -Jl  8  (section  8.2) 


°fck 

6 

* 

; 

c 

l 

t 


fixed  constant,,  defined  by  (12) 
fixed  constant,  defined  by  (16) 

geocentric  latitude  (declination);  e/(1  ♦  q)  in  section  3 
symbol  fo*’  pure  short-pe.  iod  perturbation  ,(5p  in  Part  1), 
generic  orbital  element  (osculating) 
mean  element  corresponcing  to  c 

serai -mean  element 

rate  of  change  of  due  to 
'-Or 
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K 

X 


u 


P 


0 


l 


Uk 


* 


Q 


/ 


LIST  OF  SYMBOLS  (concluded) 

Index  related  to  k  ,,  but  of  opposite  parity 
-Jl  R/p  (section  8.2) 

Earth’s  gravitational  constant 
quantity  such  that  p  -  o  ♦  qip 
modified  mean  anomaly  at  epoch 
summation  (different  use  of  I  in  Part  1) 

1  if  k  -  0,  2  if  k  >  0 
quantity  such  that  $  •  »  ♦  ch 
argument  of  perigee 
u  -  iff 

right  ascension  of  the  node 
jn  dt 
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